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Аннотация 

Введение. Математическое моделирование эффективно при анализе промышленной безопасности на 
металлургических предприятиях, в частности для отслеживания проблем системы «человек — машина». Чтобы 
ввести в рассмотрение фактор времени, задействуют рекуррентные соотношения (в дискретной модели) и 
дифференциальные (в непрерывной). Однако необходимо также решить проблему привязки параметров модели 
к реальным условиям производственной среды и к человеческому фактору. Цель данного 
исследования — создание метода определения параметров имитационных математических моделей динамики 
психофизиологических показателей оператора, влияющих на его работу. 

Материалы и методы. Психофизиологическое состояние (ПФС) оператора оценивали по работоспособности, 
утомляемости и ошибаемости. Данные собрали по тесту цифровой корректурной пробы (ЦКП). На основании 
полученных результатов вычислили экспериментальные значения показателей ПФС оператора, которые привели 
к нормированной шкале [0, 1]. Эти показатели для конкретного респондента, математическую модель и 
разработанный алгоритм задействовали при определении числовых значений параметров модели. Для 
интерпретации показателей работоспособности, утомляемости и ошибаемости ввели шкалы с пятью градациями. 
Результаты исследования. Использование модифицированного авторами варианта математической модели 
показало значительное улучшение ее прогностических свойств. Из 10 участников наилучший результат оказался у 
респондента № 7, худший — у респондента № 8. В течение 1-го часа работы (с 9.00 до 10.00) их работоспособность 
выросла примерно одинаково, с 0,5-0,55 почти до 0,6. Затем показатель респондента № 7 активно увеличивался и до 
конца рабочего дня оставался существенно выше уровня «хороший». Показатель респондента № 8 падал и с 14.00 до 
15.00 оказался ниже среднего. Разницу во многом определили хронотипы операторов. Их хронофизиологические 
особенности сказались также на утомляемости и ошибаемости. Для разных участников экспериментов 
варьировалось качество модели. В одном случае оно оказалось отличным (средняя относительная ошибка <5 %), в 
трех случаях — хорошим (<10 %), в четырех — удовлетворительным (<15 %). 

Обсуждение и заключение. Предлагаемый подход дает возможность получить для каждого индивидуума 
динамические профили его психофизиологических характеристик, оценить их взаимосвязи и выполнить прогноз 
на основе модифицированной математической модели. Однако для расширения функциональных возможностей 
моделей в реальных условиях работы оператора металлургического оборудования требуется увеличить объем 
выборки, уменьшить шаг дискретного времени и выполнить исследования для различных условий работы с 
учетом технологических, климатических, экологических, психологических и прочих факторов. 


Ключевые слова: безопасность на металлургических предприятиях, система «человек —[ машина», 
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Введение. Согласно данным Федеральной службы по экологическому, технологическому и атомному 
надзору Российской Федерации, к началу 2022 года на металлургических и коксохимических производствах 
зарегистрировано 1280 опасных производственных объектов, 19 из которых относятся к первому классу 
опасности, 325 — ко второму, и 922 — к третьему!. Как пишет журнал «Промышленные страницы»?, в 
2022 году на российских металлургических предприятиях самыми частыми происшествиями были пожары и 
взрывы. Эти инциденты обусловлены техническими причинами и нарушениями правил промышленной 
безопасности. Например, в январе 2022 года трое сотрудников Новолипецкого металлургического комбината 
проводили плановый осмотр межцехового газопровода без обязательной газозащитной аппаратуры. Рабочие 
погибли в результате отравления токсичными парами и газом. 

В [П выделены технические, санитарно-гигиенические, организационные и психофизиологические группы 
факторов аварийности и травматизма на металлургических и коксохимических предприятиях. Последние две 
группы тесно связаны с человеком — участником производственного процесса. С точки зрения 
производственной безопасности металлургических предприятий особая ответственность возлагается на 
операторов высокотехнологичных агрегатов. В инженерной психологии их принято называть 
«человек-оператор». В рамках представленной работы в этом значении будем употреблять слово «оператор». 
Иногда именно он предотвращает переход опасности в инцидент или аварию на производстве. 

ГОСТ 12.0.003-20153 выделяет в самостоятельный блок опасные и вредные производственные факторы 
психофизиологического воздействия на человека. Оператор испытывает нервно-психические перегрузки, 
связанные с напряженностью трудового процесса. В их числе: 

— умственное перенапряжение; 

— перенапряжение анализаторов; 

— монотонность труда; 

— неуверенность в действиях из-за недостатка образования и опыта. 

Эмоциональные перегрузки ведут к переутомлению, плохому самочувствию, стрессовому состоянию ит. п. [1]. 

Надежность человека как элемента сложной технической системы зависит от внутренних и внешних 
условий, которые меняются во времени. Есть еще одна переменная — сам человек. Работу выполняют люди с 
разными личностными качествами, здоровьем, опытом и пр. 

В [2] исследовались сведения о несчастных случаях (НС) в кузнечно-прессовом производстве (КПП) за 
1985, 1987 и 1989 гг. Отмечены некоторые особенности в распределении по количеству часов от начала работы 
до инцидента. Эти закономерности можно объяснить человеческим фактором, точнее суточным ритмом 
трудовой деятельности (рис. 1). 
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Рис. 1. Распределение НС в кузнечно-прессовом производстве по количеству часов от начала работы до НС [2] 


' Годовой отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору в 2021 году. Москва: 2022. 
С. 407. ОВГ: Би рз://\ууу\.созпадтог.га/рябНс/аппиа! терог5/Годовой%200отчет%20за%202021%20г.р4Е (дата обращения: 22.11.2023). 

? Опасная тенденция: аварии и ЧС на металлургических производствах в 2022-2023 годах. Промышленные страницы. ОВГ: 
Брз:/Лп4разез.га/за#е/ауагееее-па-туеаПигоееснуезКееВ-ргоеехуо4зуа/?узсПа=шз$696072208396380 (дата обращения: 22.11.2023). 

3ГОСТ 12.0.03—2015. Система стандартов безопасности труда. Опасные и вредные производственные факторы. Электронный фонд 
правовых и нормативно-технических документов. ОВГ: В рз://40сз.сп.га/4оситепИ 1200136071 (дата обращения: 21.11.2023). 
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На графиках рис. 1 достаточно выражена синхронность пиков, а также минимальная аварийность в 
последние часы работы. Значит, время работы оператора — существенный фактор возникновения несчастного 
случая на производстве. 

В [3] экспериментально установлена необходимость учитывать хронотип оператора при определении 
профессионально важных качеств. 

Для предприятий металлургической отрасли методы математического моделирования эффективны при 
анализе ситуаций, которые: 

— связаны с промышленной безопасностью; 

— могут возникнуть в сложной технической системе «человек — машина»; 

— формируются и развиваются в различных условиях внутренней и внешней среды. 

Математические модели, описывающие психофизиологическое состояние оператора, разрабатывались в 
различных математических постановках. При этом задействовали инструментарий измерения первичных 
показателей. В.Г. Абашин выяснил, как использование биометрических технологий связано с 
психофизиологическим состоянием оператора автоматического рабочего места, его работоспособностью, 
снижением количества брака и аварий технологического процесса по вине человека [4]. Работоспособность 
оператора моделируется по клавиатурному почерку. С этой целью применяют нечеткие множества как основу 
адаптивных моделей — биометрической и мультибиометрической [5]. 

Авторы [6] используют в математической модели биометрические признаки (голос, клавиатурный почерк и 
характер работы с компьютерной мышью). Это позволяет судить о психофизиологическом состоянии 
оператора: норма, усталость, опьянение, возбуждение, расслабление (засыпание). Модели базируются на 
стратегии Байеса, а также нейросетевом подходе и позволяют оценить уровень ПФС оператора и 
спрогнозировать его способность к выполнению текущих операций. 

В [7] разработаны имитационные математические модели взаимосвязи факторов психофизиологического 
состояния оператора. Авторы задействовали рекуррентные соотношения в дискретной модели и системы 
дифференциальных уравнений в непрерывной модели, что позволяет ввести в рассмотрение фактор времени. В 
этом случае параметры модели определяются по: 

— результатам тестирования конкретного работника; 

— его способности к определенному виду деятельности; 

— итогам экспериментальных исследований рабочего места и функциональных действий оператора. 

Эта часть моделирования меньше других формализуется и в наибольшей степени влияет на возможности 
использования модели для корректировки реальных действий оператора. В настоящее время проблема не 
решена. Таким образом, сохраняют актуальность вопросы определения параметров математических моделей, 
которые могут корректно описать взаимосвязи различных факторов психофизиологического состояния 
оператора. Цель данной статьи — разработка методики определения параметров таких моделей. Научные 
изыскания в данной сфере будут способствовать выявлению показателей состояния оператора 
металлургического оборудования, которые влияют на его функциональные характеристики и могут стать 
причиной серьезных ошибок с негативными последствиями. 

Материалы и методы. Для оценки психофизиологического состояния оператора выбрали следующие 
показатели: 

— работоспособность (возможность целенаправленно выполнять работу определенного количества и 
качества за обозначенное время); 

— утомляемость (снижение работоспособности с нарушением координации движений, падением 
концентрации внимания и точности решений [8]); 

— ошибаемость (оценивается количеством ошибочных действий). 

Утомление определяется по субъективным и объективным признакам [9]. Субъективный уровень утомления 
устанавливали по шкале оценки усталости“ (ее аззеззтепё зсае, РАЗ) [10]. Ее разработала группа 
нидерландских ученых под руководством Х. Михильсен (Н.Т. МасШе|5еп). Это 10 вопросов о ежедневном 
самочувствии респондента. 

В качестве источника объективных параметров в данной работе использовали тест цифровой корректурной 
пробы (ЦКП)5. По его результатам  вычислялись экспериментальные значения показателей 
психофизиологического состояния оператора. Числовые значения параметров модели определяются по 
названным показателям для конкретного респондента в совокупности с используемой моделью. 


* Шкала оценки усталости: Еабеме Аззеззтепе Зсае (ЕАЗ). ПВГ: В&рз:/Лмазое.ого/дупапис/те(#а/78/Чоситет/Опезбопайтез/ аз газ апоп.Ва (дата 
обращения: 21.11.2023). 
5 Цифровая корректурная проба. ОВТ.: № рз:/лтею4ог.га/ Че /когтеКкотгекисВ15.рНр (дата обращения: 18.08.2023). 
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Задача подбора параметров математической модели взаимосвязи различных психофизиологических 
факторов оператора решается в два этапа. 

1-й — экспериментальный. Участников тестируют, собирают данные об условиях работы оператора и его 
личностных качествах, которые могут влиять на исследуемые параметры трудовой деятельности. 

2-й — расчетный. Используется база данных, собранных на 1-м этапе. С учетом выбранной математической 
модели оцениваются ее параметры для каждого конкретного работника. При достаточном объеме 
статистических данных в широком диапазоне варьируемых параметров строятся усредненные модели, которые 
можно применить для прогнозирования изменения психофизиологических показателей работы оператора в 
различных условиях и с учетом временного фактора. 

В данном исследовании для электронного тестирования использовали программу, представленную 
интернет-ресурсом «Интерактивный портал — книга методик саморазвития и достижения успеха»°. 

При выполнении ЦКП испытуемый строка за строкой просматривал сгенерированный программой числовой 
массив и зачеркивал цифры, указанные в задании. На выполнение теста давали 3 минуты. Затем программа 
выдала описанные ниже результаты. 

1. Основные первичные показатели: 

— время выполнения теста 1; 

— общее количество просмотренных цифр до последней выбранной цифры № 

— общее количество просмотренных строк С; 

— общее количество цифр, которые необходимо вычеркнуть и; 

— общее количество вычеркнутых цифр М; 

— число верных ответов 5; 

— число пропущенных цифр Р; 

— число ошибочно выбранных цифр О. 

2. Расчетные показатели, характеризующие: 

— скорость (производительность) внимания 4; 

— точность работы 7 (в трех вариантах); 

— умственную продуктивность Ё; 

— умственную работоспособность Ан; 

— концентрацию внимания К; 

— устойчивость концентрации внимания К»; 

— объем зрительной информации И; 

— скорость переработки информации 0. 

Для решения поставленной задачи отобрали показатели, на основе которых оценивались 
работоспособность Х, утомляемость Уи ошибаемость 7 оператора по формулам (1-4). 


1. Умственная работоспособность А, (по цифровой корректурной пробе”, единица измерения — знаков в секунду): 
А, =(№М/1)х((М-(О+Р))/п). (1) 


Для приведения показателя к нормализованной шкале [0, 1] в безразмерных единицах авторы предлагают 
формулу (2). Результатом будет 1-тое значение показателя работоспособности Х' (соответствует дискретному 
условному времени 1: 


х; = Аи; [тах { А, } > (2) 


где тах{А„} — максимально возможное значение умственной работоспособности по данному виду теста 
независимо от тестируемого работника. 

Рассмотрим стандартный тест ЦКП с 1600 цифрами. Допустим, респондент за 180 секунд полностью и 
правильно выполнит тест. В этом случае величина тах{А,} составит 8,889 знаков в секунду. Если же 
тестируемый выполнит тест менее чем за 180 с, значение Х при расчете по формуле (2) может оказаться 
больше 1. Тогда для Х принимается значение 1. 

2. Утомляемость У оценивается по показателю К (концентрация внимания): 


У=(1-К/100), (3) 


где К = (М-0) п 100 — коэффициент, характеризующий концентрацию внимания, %. 


° Интерактивный портал — книга методик саморазвития и достижения успеха. ОВТ.: № рз://тею4от.ги/ (дата обращения: 18.10.2023). 
7 Цифровая корректурная проба. ОВГ.: В рз://теюогЕ ти ез 5 ЛкоггекИКоттекс15.рНр (дата обращения: 18.10.2023). 
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Данный показатель У принимает безразмерные значения в интервале [0, 1]. 
3. Ошибаемость 7 предлагается оценить по формуле: 


2=(0-+Р)/(М+Р). (4) 

Данный показатель также принимает значения в интервале [0, 1]. 
В дальнейшем будем характеризовать Х, У и 1 как унифицированные количественные индикаторы 
психофизиологического состояния оператора, связанные с его основной производственной деятельностью. Для 


трактовки показателей Х, Уи 7 разработана шкала, представленная в таблице 1. 


Таблица 1 
Интерпретация значений показателей в качественной форме 
Показатель Интервалы значений показателей / качественная характеристика 
[0-0,20] [0,20-0,40] [0,40-0,60] [0,60-0,80] [0,80-1,00] 
Работоспособность Х 
низкая ниже средней средняя хорошая высокая 
[0-0,20] [0,20-0,40] [0,40-0,60] [0,60-0,80] [0,80-—1,00] 
Утомляемость У 
низкая ниже средней средняя выше средней высокая 
[0-0,01] [0,01-0,05] [0,05—0,10] [0,10-0,20] [0,20-1,00] 
Ошибаемость 7 
незначительная заметная существенная значительная высокая 


Показатели Х, Уи 7 — безразмерные, унифицированные со значениями в интервале [0, 1]. 

На втором этапе для подбора параметров математической модели есть два способа. 

1. При помощи модуля «Поиск решения» программы Ехсе| создается сценарий решения системы для 
некоторых начальных параметров. В его основе — система рекуррентных уравнений математической модели в 
дискретной форме. Затем «Поиск решения» выдает оптимальные значения параметров модели для целевой 
функции, равной сумме квадратов отклонений расчетных значений показателей Х, У и Г от экспериментальных. 

2.В программе Апу|о21с*, которая позволяет в автоматическом режиме подбирать параметры 
рассматриваемых математических моделей в дискретной форме рекуррентных соотношений и в виде системы 
дифференциальных уравнений. Апу[021с целесообразно использовать после анализа результатов 
предварительной оценки параметров модели с помощью модуля «Поиск решения» Ехсе]. 

Для апробации предлагаемой методики отобрали 10 человек — студентов и сотрудников Донбасского 
государственного технического университета. Они работали как операторы на учебных тренажерах с 
автоматическим и полуавтоматическим управлением. В начале и в конце рабочего дня оценивался 
субъективный уровень утомления по шкале ЕРАЗ. В начале каждого часа, с 9.00 до 16.00, выполнялись тесты 
ЦКП. Данные обрабатывали по формулам (1-4) для каждого респондента. Результаты тестирования и их 
унификацию (Х, Уи 2) использовали на втором этапе для определения параметров математической модели. 


В данной работе применили модель в рекуррентной форме [7], которая описывается системой уравнений: 


Ха = № + = а т) Вх, у; — Йх,2, 
4х; + с! [4 
У; У; 
= + (1-29) +Бжу,, 5 
Ян =У 2 ру +, ( ь) 2%: У (5) 
7. РАВ 
22, +а ——.(-—^) +6: у,2,. 
+ 2. +6, ( ы) ву 


Решением системы (5) являются три сопряженных временных ряда длиной т индикаторов 
Х = {хь 1=0, 1,2,..., т}, У= {5 1=0, 1,2,....т}, ПД = {2 1=0, 1,2,...,т}. Индекс: — переменная, 
характеризующая дискретное время в системе. Параметры а» В, К, Й, 4, с, (=1,2,3) определяются в 


результате решения оптимизационной задачи по исходным данным тестирования. 


$ АпуГозтс: имитационное моделирование для бизнеса. ОВТ.: Вирз://\у\у.апуюз1с.ги/ (дата обращения: 12.11.2023). 
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Математическая модель оптимизационной задачи описывается системой ограничений (6) и целевой 
функцией (7). 
бы(арь р, ьй,Ч су) = (Р-Р У 54 
в 
би(арь р, ь в, ,с,)= (289 -2Р <, 
аб „Кв, Ас; >0, (6) 
хР > 0, у? >0,2Р > 0, 


ХР < Ту? <1,2Р <1, 


1=1,...т, 
1=1,2,5. 
(а ,ь „и, 4,с,)=У`8и+У ви + Уи тт. (7) 
1 1 = 


Здесь индексом ф помечены фактические значения индикаторов Х, Уи 7, найденные в результате обработки 
данных тестирования оператора, р — расчетные значения Х, Уи 7, определяемые из решения системы (5) и 
зависящие от искомых в данной задаче оптимизации переменных а, Б,, К, й, 4, с. (=1, 2, 3). На варьируемые 
переменные а, Б, К, й, Ц, с; ограничения задаются по максимально допустимому значению А квадратов 
отклонений (6х, бу, 6-2) расчетных значений Х, Уи 1 от фактических для всех т значений рассматриваемых 
временных рядов. Целевая функция (а, Б,, К, й, а, с}) равна сумме квадратов отклонений (бл, бук, 0). 

Для решения оптимизационной задачи применялся метод обобщенного приведенного градиента (ОПГ). 

Результаты исследования. 1-й этап — экспериментальный. Результаты тестирования участников 
эксперимента (респонденты №1-№10), преобразованные к унифицированным показателям Х, Уи 7, представлены 
на диаграммах размаха (рис. 2-4). В виде точки показан центр распределения (медиана), прямоугольник 
обозначает границы варьирования (квартили 25 % — 75 %), «усики» — нижнюю и верхнюю границу значений 


показателя (ти — тах). Диаграммы выполнены с помощью программного комплекса «ЗайзИса». 
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Рис. 2. Показатели работоспособности операторов 
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Рис. 3. Показатели утомляемости операторов 
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Рис. 4. Показатели ошибаемости операторов 


Анализ рис. 2-4 позволил выделить респондентов с лучшими и худшими оценками их деятельности 
(таблица 2). Учитывались уровни показателей Х, Уи 7, соотнесенные с градациями таблицы 1, а также квартили 


и интервалы варьирования показателей. 
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Таблица 2 
Ранжирование респондентов по результатам тестирования 
. Варианты (№ респондента) 
Критерий оценки 
наихудшие наилучшие 

Работоспособность Х (чем больше, тем лучше) 8, 1, 4,6 7, 9, 3 
Утомляемость 7 (чем больше, тем хуже) 8, 1, 4,6 7, 10 
Ошибаемость 7 (чем больше, тем хуже) 4,8, 1 9, 10,7 
Интегрально (с учетом Х, У, 7) 8,4, 1 7, 10, 9 


Графики на рис. 5 отражают динамику работоспособности для респондентов №7 (пунктирная линия, 
лучший результат) и № 8 (сплошная линия, худший результат). 
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Рис. 5. Динамика работоспособности операторов 


Для объяснения достаточно большой разницы данных респондентов 7 и 8 учитывали их возраст, пол, род 
занятий, самочувствие и др. Наиболее существенным оказался фактор хронотипа. Классификация, принятая в 
1970 году, подразумевает три категории людей с отличительными поведенческими признаками и генетической 
разницей биоритмов. Это так называемые «жаворонки», «совы» и «голуби» [11]. Респондент № 8 по 
хронотипу — «жаворонок», №7 — «сова». В результате работоспособность № 8 снижается к концу дня, 
а №7 — повышается (рис. 5). Хронофизиологические особенности операторов проявляются также в динамике 


их утомляемости и ошибаемости. 


2-й этап — расчетный. Получены параметры модели (5), рассчитанные для каждого участника 
эксперимента. С этой целью решили оптимизационную задачу (6)-(7). Диапазоны изменения параметров 


приведены в таблице 3. 


Таблица 3 


Диапазоны параметров математической модели (5) для группы участников эксперимента 


Показатель Диапазоны параметров Показатель Диапазоны параметров 

а 0,0004—0,0300 Бз 0,0001—10,1400 
Ь 0-—0,9700 Аз 0,0020-0,0700 
К 0,5900-—1,0000 С 0,0010-0,2000 
р 0—1,8100 © 0,0001-0,0010 
а2 0-—0,0010 сз 0,0100-0,3300 
Ь> 0,0070-0,2600 4 =0 

2 0,4200—1,0000 р =0 

аз 0,0005-—0,0600 4 0,0400-0,9800 


В таблице 4 представлены результаты оценки качества математических моделей динамики изменения 
показателей Х, Уи 2, оцениваемых по среднеквадратической ошибке модели и средней относительной ошибке 


для всех испытуемых. 
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Таблица 4 
Оценка качества математических моделей 


№ респондента 1 2 3 4 5 6 7 8 9 10 


Средняя относительная ошибка =, % | 14,27 | 4,54 | 20,93 | 12,80 | 7,28 | 12,00 | 9,66 | 15,03 | 12,23 | 9,70 


Среднеквадратическая ошибка, ед. | 0,182 | 0,107 | 0,253 | 0,267 | 0,158 | 0,247 | 0,271 | 0,256 | 0,192 | 0,166 


Для разных респондентов показатели качества модели варьируются в широком диапазоне. Отличное 
качество (&<5%) получено в одном случае, хорошее (=<10%) — в трех случаях, удовлетворительное 
(=<15%) — в четырех. Для двух случаев (№3 и №8) не удалось удовлетворительно решить задачу 
оптимального подбора параметров математической модели. Это говорит о том, что не все влияющие факторы 
учтены либо модель (5) в некоторых случаях не работает. 

Согласно результатам анализа алгоритма решения поставленной задачи, ряд упрощений улучшит 
сходимость результатов. Можно, например: 

— уменьшить на 3 ед. число параметров системы (5), сократив числители и знаменатели вторых слагаемых 
правых частей системы на ст, со, сз соответственно; 

— удалить из системы (6) ее первое, второе и третье ограничения, переведя их в статус наблюдаемых ограничений. 

В итоге алгоритм метода ОПГ будет работать лучше. 

На рис. 6 и 7 приведены сравнительные диаграммы оценки качества построенных математических моделей 
в двух вариантах: 

—1 — исходная модель (5); 

— 2 — модифицированная модель. 

Как видно из диаграмм, решение задачи подбора параметров математической модели динамики показателей 
ПФС оператора значительно улучшено при использовании модели 2 для всех случаев, кроме № 8. 
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Рис. 6. Сравнение математических моделей по средней квадратической ошибке 
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Рис. 7. Сравнение математических моделей по средней относительной ошибке 
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Обсуждение и заключение. Предложенный подход позволяет получить и прогнозировать для каждого 
индивидуума динамические профили его психофизиологических характеристик, составить математические 
модели взаимосвязей. Целесообразно увеличить объем выборки, уменьшить шаг дискретного времени и 
выполнить исследования для различных условий работы, включающие технологические, климатические, 
экологические, психологические и прочие факторы, которые могут влиять на эффективность работника. Это 
следует сделать для расширения функциональных возможностей моделей, их применения в реальных условиях 
работы оператора металлургического оборудования. 

Обнаруженный эффект влияния хронотипа человека на его ПФС требует более основательных 
исследований. В статье [12], посвященной поиску компонентов циркадных часов у людей, на больших 
статистических выборках показано, что хронотип зависит от многих факторов (пол, возраст, график работы и 
др.). Все это необходимо диагностировать и учитывать в математических моделях оценки ПФС операторов. 

Интересно также дополнить модель компонентами культуры производственной безопасности [13], которые 
можно привести к количественному виду и использовать как поправочные коэффициенты. 

В дальнейшем предполагается усовершенствовать модель дополнительными переменными. Для этого нужно 
зафиксировать психофизиологические показатели и определять местонахождения оператора в режиме 
реального времени. 

Представленную в работе методику можно взять за основу при решении описанных ниже задач. 

— Составление и анализ динамического профиля работника при приеме на работу в качестве оператора 
машины, агрегата или устройства, где важны такие характеристики ПФС, как работоспособность, утомляемость 
и ошибаемость. Эта задача решается для конкретного человека в определенных условиях производства, что 
позволяет подобрать оптимальный режим работы и отдыха оператора, способствует сохранению здоровья 
специалиста и повышению уровня промышленной безопасности предприятия. 

— Формирование и ведение специализированных баз статистических данных, включающих характеристики 
рабочих мест и сотрудников, в том числе их психофизиологическое состояние. На основе собранного 
статистического материала и систем мониторинга ПФС оператора в режиме реального времени можно строить 
прогнозные модели для предупреждения нештатных и аварийных ситуаций на металлургических предприятиях. 
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Аннотация 

Введение. Нормативными документами по пожарной безопасности при определении численности и 
технической оснащенности подразделений пожарной охраны в явном виде не учитываются природно- 
климатические и географические особенности субъектов Российской Федерации. При этом в научной 
литературе встречаются сведения о влиянии отдельных природно-климатических факторов на эффективность 
действий пожарной охраны при тушении пожаров. Целью настоящего исследования явилось определение 
влияния системы природно-климатических и географических характеристик субъектов Российской Федерации 
на показатели оперативного реагирования территориальных подразделений пожарной охраны. Полученные 
результаты рекомендованы для дальнейшего использования при нормировании численности и технической 
оснащенности подразделений пожарной охраны. 

Методы и материалы. Проведен факторный анализ статистических данных за 2020—2022 гг. Статистические 
данные получены из федеральной государственной информационной системы «Федеральный банк данных 
«Пожары», с сайта Федеральной службы государственной статистики Российской Федерации и из других 
источников. Для анализа отобраны 10 показателей, характеризующих природно-климатические и 
географические особенности субъектов Российской Федерации, 2 показателя оперативного реагирования 
территориальных подразделений пожарной охраны и 4 показателя обстановки с пожарами. 

Результаты исследования. Выделены пять значимых факторов, изменение которых объясняет изменение 
наблюдаемых показателей. Первый фактор характеризует связь климатических условий с показателями 
обстановки с пожарами. Второй фактор связывает показатели оперативного реагирования подразделений 
пожарной охраны с особенностями рельефа субъектов Российской Федерации. Третий фактор описывает 
взаимосвязь показателей обстановки с пожарами и показателей оперативного реагирования с плотностью 
населения и лесистостью территории. Остальные факторы в показатели обстановки с пожарами и оперативного 
реагирования существенного вклада не вносят. 

Обсуждение и заключение. С помощью математического анализа и с использованием факторного 
моделирования авторами исследована взаимозависимость природно-климатических и географических 
особенностей субъектов Российской Федерации, показателей обстановки с пожарами и показателей 
оперативного реагирования подразделений пожарной охраны. Определены наиболее значимые факторы, 
влияющие на эти показатели. Среди них средняя температура воздуха, лесистость территории, наличие горных 
массивов, плотность населения. Данные показатели следует учитывать при определении численности и 
технической оснащенности подразделений пожарной охраны для повышения эффективности их 
функционирования. 


Ключевые слова: пожарная охрана, факторный анализ, субъект Российской Федерации, гибель, травматизм, 
время прибытия, природно-климатические и географические показатели 
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Введение. При определении численности и технической оснащенности территориальных подразделений 
пожарной охраны проводится обследование территории субъекта Российской Федерации, которое включает в 
себя анализ состояния существующей системы обеспечения пожарной безопасности субъекта Российской 
Федерации, оперативно-тактических характеристик (особенностей) рассматриваемой территории, частоты 
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возникновения пожаров, рисков гибели и травмирования на них людей'. Среди факторов, влияющих на 
состояние пожарной безопасности в субъектах Российской Федерации, большое значение имеют природно- 
климатические и географические характеристики. В этой связи в настоящей работе с использованием 
факторного анализа было проведено исследование взаимосвязи природно-климатических и географических 
характеристик субъектов Российской Федерации, показателей оперативного реагирования территориальных 
подразделений пожарной охраны и показателей обстановки с пожарами. Факторный анализ является одним из 
наиболее широко используемых методов, направленных на выделение из множества разнородных показателей 
небольшого количества латентных факторов, вариация которых может объяснить изменение большей части 
наблюдаемых показателей [1-5]. Это позволяет построить математическую модель, в которой исследуемые 
факторы имеют простую и наглядную структуру [6]. Ранее авторами с использованием факторного анализа 
проводилось изучение оперативного реагирования подразделений пожарной охраны [7], социальных 
последствий пожаров [8], обстановки с пожарами в сельских поселениях [9], готовности подразделений 
пожарной охраны [10]. 

Материалы и методы. Проводимое исследование включало в себя три этапа. На первом этапе была 
сформирована матрица исходных показателей. Для этого в каждом субъекте Российской Федерации были 
определены значения 16 показателей, разбитых на три группы. 

В первую группу вошли данные, характеризующие обстановку с пожарами (обозначение показателя указано 
в скобках): 

— количество пожаров, приходящихся на 1 тыс. человек населения, ед. (71); 

— средний прямой материальный ущерб от одного пожара, руб. (75); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Уз); 

— количество травмированных при пожарах, приходящихся на 100 тыс. человек населения, чел. (74). 

Во второй группе отражались показатели оперативного реагирования подразделений пожарной охраны: 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— среднее время, затрачиваемое на тушение пожара, мин. (22). 

Третью группу составляли показатели, характеризующие природно-климатические, географические и 
демографические особенности субъектов Российской Федерации: 

— плотность населения, чел/км? (А); 

— доля водной поверхности от площади территории субъекта, % (>); 

— длина морской береговой линии от общей протяженности границы субъекта, % (ХЗ); 

— доля лесов от площади территории субъекта, % (Ха); 

— среднегодовое количество осадков, мм в год (Х°5); 

— средняя температура июля, °С (Х5); 

— средняя температура января, °С (Х7); 

— сейсмическая опасность (наличие населенных пунктов с указанной сейсмической интенсивностью для 
степени сейсмической опасности), °С (Хз); 

— доля территории, занятой горами, от общей площади территории субъекта, % (Хэ); 

— количество солнечных дней в году, (По). 

Показатели обстановки с пожарами и показатели оперативного реагирования подразделений пожарной охраны 
получены из электронной базы данных пожаров, которая ведется с использованием автоматизированной 
аналитической системы поддержки и управления контрольно-надзорными органами МЧС России?. 

Сведения о плотности населения и площади лесов в субъектах Российской Федерации взяты из данных 
Федеральной службы государственной статистики Российской Федерации?. Характеристики климата субъектов 
Российской Федерации получены на справочно-информационном портале «Погода и климат»“. При 
определении сейсмической опасности в субъектах Российской Федерации использованы данные 
СП 14.13330.20145. 


1 Об утверждении методик расчета численности и технической оснащенности подразделений пожарной охраны. Приказ МЧС России от 
15.10.2021 № 700. ОВГ: В#рз://базе.сагапй.ги/403136953/ (дата обращения: 18.12.2023) 

206 утверждении Регламента работы в информационной системе «Автоматизированная аналитическая система поддержки и 
‘управления контрольно-надзорными органами МЧС России». Приказ МЧС России № 954 от 04.10.2022. ОВГ: Врз://Ягетап.с1а6/погтаНнуе- 
Ч4оситеп/риКа7-псВ$-гоз31-954-01-04-10-2022-оБ-шуегив4епи-гезатеша-габобу-у-шРппас71оппо]-515ете/ (дата обращения: 11.12.2023) 

3 Регионы России. Социально-экономические показатели 2022. Федеральная служба государственной статистики Российской Федерации. 
ОВГ: Би рз://тозз{а. соу ги/КоЧет/2 10/4осиеп 13204 (дата обращения: 08.12.2023) 

* Справочно-информационный портал «Погода и климат». ОВГ.: ВИр://\у\у\ .росодай таги (дата обращения: 08.12.2023) 

5 Свод правил. Строительство в сейсмических районах. СП 14.13330.2018 от 25.11.2023. ОВГ.: № рз://Аосз.сп9.гл/4оситеп/ 550565571 (дата 
обращения: 08.12.2023) 
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Отобранные статистические показатели характеризовались болыпим разбросом, что было обусловлено 
различными размерами субъектов Российской Федерации (площадь территорий и численность населения), а 
также различиями их природно-климатических и географических условий. В связи с этим для приведения 
статистических данных к единообразной форме была проведена нормировка ряда показателей, что позволило 
преобразовать множество натуральных показателей в синтетические. Нормировка показателей осуществлялась 
на площадь территории субъектов (плотность населения, доля лесов и водной поверхности от площади 
территории субъекта) и на численность населения (количество пожаров, число погибших и травмированных 
при пожарах). 

Дальнейшее моделирование проводилось с использованием полученной матрицы синтетических 
показателей. Значимые факторы выделялись с использованием трех наиболее распространенных методов: 
центроидного, главных компонент и максимального правдоподобия. Все три метода дают близкие результаты, 
которые в пределах статистических погрешностей согласуются между собой. В дальнейшем для проведения 
факторного анализа был использован метод главных компонент. 

На втором этапе с помощью метода главных компонент были получены значимые факторы. Для облегчения 
их интерпретации выполнено вращение полученных факторов с использованием следующих методов: 

— варимакс (данный метод путем уменьшения числа переменных для каждого фактора позволяет выполнить 
лучшее разделение факторов); 

— квартимакс (данный метод путем уменьшения числа факторов, относящихся к каждой переменной, дает 
возможность выделить генеральный фактор и упростить интерпретацию); 

— биквартимакс; 

— эквимакс. 

Перед проведением процедуры вращения факторные нагрузки были нормализованы с использованием 
метода Кайзера. Это позволило исключить влияние переменных с большой схожестью на итоговый результат. 
В последующих исследованиях использовался метод вращения варимакс, поскольку анализ полученных 
результатов показал, что факторные нагрузки несущественно зависят от метода вращения. 

В результате факторного анализа было получено пять значимых факторов. Доля общей дисперсии, 
объясняемой данными факторами, распределена следующим образом: первый фактор — 22,1%, второй 
фактор — 15,4 %, третий — 12,1 %, четвертый — 9,7 %, пятый фактор — 7,7 %. Данные пять значимых 
факторов объясняют 67,0 % общей дисперсии. 

В таблице] представлены значения факторных нагрузок для каждого из пяти факторов (значимые 
показатели выделены полужирным шрифтом). 


Таблица 1 
Матрица факторных нагрузок 
Переменная Фактор 
(показатель) 1 2 3 4 5 

У 0,019 -0,179 0,681 0,212 -0,308 

У 0,692 0,028 -0,241 0,027 0,123 

Уз 0,692 0,164 0,484 -0,123 0,021 

7 0,821 -0,095 0,229 0,046 -0,121 

7 -0,362 0,587 0,417 0,091 0,306 

22 0,009 0,218 0,509 -0,343 0,265 

М -0,093 0,069 -0,587 -0,086 -0,063 

№ 0,018 0,234 0,056 0,070 -0,836 

№ 0,000 -0,373 -0,155 -0,422 -0,532 

№ 0,236 -0,305 0,536 -0,357 0,002 

№ -0,334 0,026 0,008 -0,758 0,016 

Х -0,636 0,032 -0,062 0,630 0,200 

А -0,646 0,357 -0,343 -0,144 0,040 

№ -0,073 -0,833 0,097 0,174 0,133 

№ 0,058 -0,841 0,180 -0,019 0,051 

Хо -0,437 -0,322 0,191 0,630 -0,081 
На третьем этапе была решена творческая задача, которая выходит за пределы формального 
метода — проведена содержательная интерпретация полученных факторов с использованием предметных 

терминов. 
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Результаты исследования. Установлено, что в первый значимый фактор наибольший вклад вносят 
следующие переменные: 

— средний прямой материальный ущерб от одного пожара, руб. (75); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Уз); 

— количество травмированных, приходящихся на 100 тыс. человек населения, чел. (У4); 

— средняя температура июля, °С (Х5); 

— средняя температура января, °С (Х7); 

— количество солнечных дней в году (А10). 

Среди показателей обстановки с пожарами в первый фактор наибольший вклад вносит показатель 
«количество травмированных при пожарах, приходящихся на 100 тыс. человек населения». Данный фактор 
характеризует связь климатических условий с показателями обстановки с пожарами. В тех субъектах 
Российской Федерации, где более суровые климатические условия (низкие температуры, малое количество 
солнечных дней), отмечаются более тяжелые последствия пожаров (больше погибших и травмированных, выше 
материальный ущерб). 

Данная взаимосвязь демонстрируется на рис. | и 2, где показаны зависимости количества пострадавших 
(погибших и травмированных) при пожарах, приходящихся на 100 тыс. человек населения, от средней 
температуры января. Видно, что с уменьшением средней температуры увеличивается количество пострадавших 
при пожарах. 
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Рис. 1. Зависимость количества погибших при пожарах, приходящихся на 100 тыс. человек населения, от средней 
температуры января в субъектах Российской Федерации. 
Прямая — результат аппроксимации линейной функцией у = -0,166х -+ 4,338 
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Рис. 2. Зависимость количества травмированных при пожарах, приходящихся 
на 100 тыс. человек населения, от средней температуры января в субъектах 
Российской Федерации. Прямая — результат аппроксимации линейной функцией 
у= -0,248х + 3,565 
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Для второго значимого фактора наиболее существенными являются следующие переменные: 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— сейсмическая опасность (наличие населенных пунктов с указанной сейсмической интенсивностью для 
степени сейсмической опасности С) (Хз); 

— доля территории, занятой горами, от общей площади территории субъекта, % (Хо). 

Второй фактор характеризует взаимосвязь показателей оперативного реагирования подразделений 
пожарной охраны с географическими особенностями субъектов Российской Федерации. Наблюдается важная 
зависимость — чем большую территорию субъекта занимают горы, тем меньше среднее время, за которое 
первое подразделение пожарной охраны прибывает к месту вызова. 

Распределение факторных нагрузок для трех групп показателей в плоскости факторов 1 и 2 показано на рис. 3. 
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Рис. 3. Значения факторных нагрузок в плоскости факторов 1 и 2. Зеленым цветом обозначены показатели первой группы, 
красным — второй группы, синим — третьей группы 


В третий значимый фактор наибольший вклад вносят следующие переменные: 

— количество пожаров, приходящихся на 1 тыс. человек населения, ед. (71); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Уз); 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— среднее время, затрачиваемое на тушение пожара, мин. (72); 

— плотность населения, чел/км? (Х'); 

— доля лесов от площади территории субъекта, % (Х4). 

Третий фактор описывает взаимосвязь показателей обстановки с пожарами и показателей оперативного 
реагирования с плотностью населения и лесистостью территории. Как показали ранее проведенные 
исследования, с ростом плотности уменьшается средняя площадь области обслуживания одного подразделения 
пожарной охраны, результатом чего является снижение времени прибытия первого подразделения пожарной 
охраны. С другой стороны, чем ниже плотность, тем болыше площадь территории, на которой проживает 
одинаковое количество людей, соответственно, на большей площади происходит большее количество пожаров. 
Интересно также отметить, что с увеличением плотности населения снижается количество пострадавших при 
пожарах, приходящихся на 100 тыс. населения. 

На обстановку с пожарами также влияет площадь лесов. С увеличение доли территории субъекта 
Российской Федерации, занятой лесами, увеличивается количество пожаров, приходящихся в расчете на 
1 тыс. человек. 

Получена зависимость количества пожаров, приходящихся на 1 тыс. человек, от плотности населения. Как 
видно на рис.4, количество пожаров уменьшается с ростом плотности населения. Данная зависимость 
описывается экспоненциальной функцией. 
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На рис. 5 представлено распределение факторных нагрузок для трех групп показателей в плоскости 
факторов 1 и 3. 
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Рис. 5. Значения факторных нагрузок в плоскости факторов 1 и 3. Зеленым цветом отмечены показатели первой группы, 
красным — второй группы, синим — третьей группы 


Для четвертого значимого фактора наиболее важным являются следующие переменные: 

— среднегодовое количество осадков, мм в год (Х5); 

— средняя температура июля, С° (№5); 

— количество солнечных дней в году (Хо). 

Данный фактор связан с климатическими характеристиками субъектов Российской Федерации. Показатели 
первой и второй групп значительного вклада в данных фактор не вносят. Следует отметить, что взаимосвязь 
климатических характеристик субъектов с показателями обстановки с пожарами хорошо прослеживается из 
анализа факторных нагрузок для первого фактора (см. выше). 


Стрельцов О.В. и др. Влияние природно-климатических и географических особенностей субъектов Российской Федерации ... 


Для пятого значимого фактора определяющими являются следующие переменные: 

— доля площади водной поверхности от площади территории субъекта, % (>); 

— длина морской береговой линии от общей протяженности границы субъекта, % (Хз). 

Данные гидрографические характеристики не оказывают существенного влияния на обстановку с пожарами и 
на показатели оперативного реагирования подразделений пожарной охраны в субъектах Российской Федерации. 

Обсуждения и заключение. Впервые с использованием факторного анализа проведено исследование 
влияния природно-климатических и географических особенностей субъектов Российской Федерации на 
уровень пожарной безопасности. Проведенный анализ показал, что наиболее существенно на обстановку с 
пожарами и деятельность пожарной охраны в субъектах Российской Федерации влияют следующие природно- 
климатические, географические и демографические показатели: средняя температура, лесистость территории, 
наличие горных массивов, плотность населения. Таким образом, при определении численности и технической 
оснащенности подразделений пожарной охраны необходимо учитывать данные особенности субъектов 
Российской Федерации. 
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Аннотация 

Введение. Внутрилабораторные сличительные испытания являются важным и актуальным мероприятием для 
обеспечения качества и достоверности результатов исследований в лабораториях. Они позволяют оценить 
точность и воспроизводимость применяемых в лаборатории методов, а также выявить возможные источники 
ошибок и несоответствий в ее работе. Результаты внутрилабораторных сличительных испытаний 
предоставляются экспертам при подтверждении компетентности в рамках аккредитации. Как правило, 
сличительные испытания проводят в лабораторных условиях в привычной и спокойной для испытателей 
обстановке. Однако лаборатории, проводящие исследования в рамках специальной оценки условий 
труда (СОУТ), вынуждены проводить выездные сличительные испытания на реальных объектах, где за 
испытаниями наблюдают заказчики, которые невольно отвлекают работников лаборатории, что напрямую 
влияет на качество измерений. Целью данной работы является оценка качества проведения выездных 
внутрилабораторных сличительных испытаний на примере определения запыленности рабочего места 
каменщика на строительной площадке и определение минимально необходимого и достаточного для этого 
количества измерений. 

Материалы и методы. Для определения запыленности рабочего места каменщика использовался весовой 
метод, который заключается в сборе пыли на фильтры с последующим взвешиванием и определением 
концентрации. Оценка качества проведения внутрилабораторных испытаний производилась согласно 
ГОСТРИСО 5725-1-2002 «Точность (правильность и прецизионность) методов и результатов измерений. Часть 
1. Основные положения и определения». 

Результаты исследования. На рабочем месте каменщика выявлено превышение разовой предельно 
допустимой концентрации пыли в 1,6 раза. Средняя концентрация пыли на исследуемом рабочем месте 
составила: К» = 9,57+0,81 мг/м3, сходимость полученных результатов, 7 = 8,68 %, относительная погрешность, 
б= 8,50 %. Выявлено, что максимально допустимая разница результатов двух испытаний составляет 0,84 мг/мз. 
Разница прямых измерений массы двух образцов должна быть не более 0,1 мг. 

Обсуждение и заключение. Полученные результаты показали возможность проведения минимального 
количества измерений, что по условиям воспроизводимости испытания признаются удовлетворительными и 
могут предоставляться экспертам для подтверждения компетентности лаборатории. В рамках СОУТ 


работникам исследуемого объекта рекомендуется использовать средства индивидуальной защиты органов 
дыхания, предохраняющие от высокодисперсной ПЫЛИ. 


Ключевые слова: сходимость, внутрилабораторный контроль, запыленность, строительная площадка, рабочее 
место, специальная оценка условий труда (СОУТ) 
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Введение. Измерения на запыленность рабочих мест проводятся в рамках специальной оценки условий 
труда (СОУТ) на строительной площадке с целью выявления и оценки вредных и опасных производственных 
факторов, а также разработки мероприятий по их устранению или минимизации [1]. Результаты специальной 
оценки условий труда являются основой для разработки и внедрения соответствующих мер по улучшению 
условий труда [2]. Это может быть изменение рабочих процессов, применение новых технологий, обновление 
оборудования или улучшение санитарных условий [3]. Важно, чтобы эти меры были реализованы с учетом 
мнения сотрудников и их потребностей [4]. 
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Проведение специальной оценки условий труда на строительной площадке имеет свои особенности, 
связанные с характером и условиями труда [5]. К таким особенностям можно отнести: 

— работу на открытом воздухе, что может привести к воздействию неблагоприятных погодных условий; 

— использование строительной техники и оборудования, которые могут создавать шум, вибрацию, пыль и 
другие загрязнения; 

— выполнение работ на высоте, что требует специального обучения и аттестации персонала; 

— работа с различными материалами и веществами, которые могут быть опасными для здоровья работников; 

— необходимость соблюдения требований охраны труда и техники безопасности, включая использование 
средств индивидуальной защиты. 

В результате выполнения СОУТ устанавливаются классы условий труда для каждой профессии, 
рекомендуются мероприятия по минимизации выявленных рисков и последствий влияния вредных 
факторов [6]. По данным СОУТ руководство организации разрабатывает мероприятия по улучшению условий и 
охраны труда, также определяются размеры дополнительных выплат для работников, занятых на производстве 
с неблагоприятными или рискованными условиями [7]. При инструментальном исследовании рабочих мест 
осуществляются замеры и анализ шума, степени запыленности, естественной и искусственной освещенности, 
температуры, влажности и пр. [8]. 

СОУТ имеют право проводить организации, обладающие соответствующей аккредитацией и опытом работы 
в данной области. Ими могут быть как специализированные организации, так и лаборатории, занимающиеся 
исследованием условий труда. 

В соответствии с Федеральным законом от 28.12.2013 №426 «О специальной оценке условий труда» 
организация, занимающаяся СОУТ, должна: 

— иметь не менее 5 штатных сотрудников с действующим аттестатом на право проведения СОУТ; 

— быть зарегистрированной в реестре организаций, проводящих СОУТ, который ведет Минтруда и 
соцзащиты РФ; 

— обладать системой менеджмента качества (СМК), соответствующей ГОСТ ГТЗОЛЕС 17025-2019 «Общие 
требования к компетентности испытательных и калибровочных лабораторий» и ГОСТР 54934-2012/ 
ОН$А$ 18001:2007 «Системы менеджмента безопасности труда и охраны здоровья. Требования»; 

— располагать аккредитованной испытательной лабораторией, которая выполняет исследования и измерения 
опасных и вредных факторов. 

Испытательная лаборатория (ИЛ) должна иметь действующий аттестат аккредитации, соответствующую 
область аккредитации, отражающую возможность исследования вредных и опасных факторов производства и 
условий труда, а также аттестованное испытательное оборудование (ИО) и поверенные средства 
измерений (СИ) [9]. Согласно постановлению Правительства РФ от 26 ноября 2021 года № 2050 «Об 
утверждении правил осуществления аккредитации в национальной системе аккредитации...» аккредитованные 
лаборатории в рамках подтверждения компетентности должны проводить сравнительные испытания для 
определения сходимости результатов измерений. Сходимость результатов измерений при проведении 
специальной оценки условий труда является одним из ключевых показателей, определяющих качество 
проводимых измерений. 

Проведенные ранее исследования в данной области затрагивали в основном межлабораторные 
сравнительные испытания (МСИ) и не рассматривали внутрилабораторные сравнительные (сличительные) 
испытания. Так, в работе Козлова А.В. приведены результаты МСИ геосинтетического материала [10]. Автор 
указывает на целесообразность нормирования требований показателей прецизионности испытаний. В 
исследовании Лекомцевой М.М. и Шендалевой Е.В. проведен анализ результатов МСИ нефтепродуктов, по 
результатам которого рекомендуется преобразовывать данные в линейную зависимость для получения 
неизменного значения стандартного отклонения [11]. В работе Ясина Дургута отмечается, что 
межлабораторные сравнительные испытания, проводимые в соответствии с ГОСТ ТЗОЛЕС 17025-2019 «Общие 
требования к компетентности испытательных и калибровочных лабораторий», выполняются более эффективно, 
по сравнению со сравнительными измерениями, организованными в соответствии с ГОСТ ТЗОЛЕС 17043 
«Оценка соответствия. Основные требования к проведению проверки квалификации» [12]. 

Между тем, согласно ГОСТ 1ЗОЛЕС 17025-2019 «Общие требования к компетентности испытательных и 
калибровочных лабораторий» весьма важны еще и результаты внутрилабораторных сравнительных испытаний, 
так как они используются для проведения корректирующих мероприятий, улучшения качества работы 
лаборатории и повышения точности результатов [13]. Поэтому целью данной работы является оценка качества 
проведения выездных внутрилабораторных сличительных испытаний. Такая оценка представлена на примере 
определения запыленности рабочего места каменщика на строительной площадке. Лаборатории, 
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осуществляющие свою деятельность в рамках СОУТ, вынуждены проводить внутрилабораторные 
сравнительные испытания на выездной основе, и условия их значительно отличаются от лабораторных. 
Сотрудники лаборатории ограничены временем допуска на объект, поэтому им важно выполнить исследования 
не только качественно, но и с минимально количеством испытаний. В связи с этим задача авторов 
статьи — определить это минимально необходимое и достаточное количество таких испытаний. 

Методы исследований. Общие санитарно-гигиенические требования к воздуху рабочей зоны, а также 
методы контроля показателей микроклимата и содержания вредных веществ в воздухе регламентирует 
ГОСТ 12.1.005-88 «Межгосударственный стандарт. (Система стандартов безопасности труда. Общие 
санитарно-гигиенические требования к воздуху рабочей зоны». Для определения запыленности использовался 
весовой метод', который заключается в сборе пыли на фильтры с последующим взвешиванием и определением 
концентрации. 

Массовую концентрацию К„ пыли определяют по формуле: 


(т„ —то)-1000 


Кн = т = , (1) 
20 
где К, —Щ_ концентрация пыли в воздухе, мг/м3; То — масса чистого фильтра, мг, т, — масса фильтра [и 
осевшими частицами пыли, мг; Го === объем воздуха, приведенный к стандартным условиям, ма. 
7, .293.Р 
= (2) 


(273+ Г)-101,33” 


где Г, — объем воздуха, прошедший через фильтр, дм3; Р — атмосферное давление, кПа; Т — температура 
воздуха на рабочем месте, °С. 

Согласно рекомендациям государственной системы обеспечения единства измерений? для обработки 
результатов определения запыленности рабочего места можно использовать метод приведения, который 
предполагает наличие значений измеряемых аргументов, полученных в результате многократных измерений. 

Результат косвенного измерения находят по формуле: 


ы А ы 
А= ) —. (3) 
= Г 
где Г — число частных значений измеряемой величины; А;_— /-е значение измеряемой величины; / — номер 


последовательного измерения от 1 до Г. 


Среднее квадратическое отклонение 5(4) случайных погрешностей результата косвенного измерения 


рассчитывают по следующей формуле: 


в (4-4) 
Уз Г(Е-1) е 
При нормальном распределении отдельных значений измеряемой величины доверительные границы 
случайных погрешностей определяют по формуле: 
А=1,.5(4), (5) 
где /, — коэффициент Стьюдента, зависящий от доверительной вероятности и числа результатов наблюдений. 


Прецизионность оценивается по показателю повторяемости (относительной среднеквадратического 
отклонения повторяемости), О, И ПО ПОКазателю точности (границ относительной погрешности), 0: 


= (6) 
где х; — 1-Й результат испытаний, полученный в условиях повторяемости; х — среднеарифметическое 
значение из и результатов испытаний в условиях повторяемости (сходимости). 

А 0, 0, 
5=—100% < 25%. (1) 
х 


' Методика измерений массовой концентрации пыли гравиметрическим методом для целей специальной оценки условий труда. МИ 
АПФД-18.01.2018. ОВГ: В рз://поппайу.Копеиг.га/доситет тоде!=1 &4оситеп 4=390372 (дата обращения: 25.11.2023). 

2 Рекомендация. Государственная система обеспечения единства измерений. Измерения косвенные. Определение результатов измерений и 
оценивание их погрешностей. МИ 2083-90. Москва: Комитет стандартизации и метрологии СССР; 1991. 
ОВГ:Ьйрз://7пауюуаг.ги/со$/2/МТ 208390 СТ Гитегешуа_Козуе.В 1 (дата обращения: 25.11.2023). 
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В процентном отношении сходимость определяют в соответствии с выражением: 


О. 
п=—^ 100 %. (8) 
х 
Для оценки сходимости результатов двух измерений необходимо использовать формулу: 
м -%| ооо 
ЕР “00%, (9) 
Хх 


где 1, х2 — результаты двух параллельных измерений. 
Результаты измерений, полученные В условиях сходимости, признаются удовлетворительными, если 
полученная СХОДИМоСсТЬ Г будет меньше или равна нормативной сходимости Гнорм: 
И, =15 %. (10) 


норм 
За результат принимается выражение (8): 
хХЕХ-А. (11) 

Минимальное количество наблюдений зависит от коэффициента вариации, т.е. отношения 
среднеквадратического отклонения к среднему значению, и выбирается предварительно из опытных данных. 
Ошибку в определении среднего значения (статистическая погрешность измерений) необходимо искать в 
предположении, что закон распределения генеральной совокупности измерений — нормальный. Зная ошибку и 
доверительную вероятность, можно однозначно установить доверительный интервал для среднего значения и 
оценить максимальную статистическую погрешность измерений как отношение половины разности границ 
ошибки среднего значения к среднему значению из опытных данных [14]. 

Результаты исследований. Проведены измерения на определение запыленности строительной площадки 
при строительстве многоквартирного жилого дома. Отбор проб производился на рабочем месте каменщика в 
зоне кладочной стены фильтрами типа АФА (аналитические аэрозольные фильтры) с помощью двух 
аспираторов ПУ-4Э, обеспечивающих отбор проб с заданным объемным расходом через поглотитель по 
четырем параллельным каналам. Для установления статистической погрешности измерений инструментами и 
сотрудниками лаборатории при заданной доверительной вероятности Р = 0,95 выполнено восемь параллельных 
измерений, результаты которых приведены в таблице 1. 


Таблица 1 
Результаты восьми параллельных испытаний 
з 55° < 
ы : = -: Ф 5 < 
ы ы ы м ы з к т о |) — 
= = = | | 
Е. & Е Е Е во еб Е : о 
Е - : Е | ЕЕ ЕЁ | 
о з я я Е Е - ЕВЕ чо 8 < Е 
5- Е Ее ^ Я 2. И 6 Я = | > 
о > т я Ф 
> |= я [6%] [) Я го © я Ф = во а © 2. 
[2 я я о о я вв я Я о = я |= 
- ПВ - а | В ЕЕ 391 Е 
2 г : Е В ЕЕ ЕВ 
о о я = я в Е 2. = В Ё ы Е 
з з 8 о Я > о 5 Е 
> > 5. = я = = = с. 5 я 
Е ЕЯ 
< ©) ©) 
п=8 
№1 30965,6 30966,5 2 102,2 60 21,1 7,46 
№2 31346,8 31347,8 2 102,2 60 21,1 8,29 
№3 33980,5 33981,6 2 102,2 60 21,1 9,12 
№4 35572,3 35573,6 2 102,2 60 21,1 10,78 | 9,33 1,39 10,29 | 14,85 
№5 34926,9 34928,3 2 102,2 60 21,1 11,61 
№6 30477,7 30479,0 2 102,2 60 21,1 10,78 
№7 32678,4 32679,5 2 102,2 60 21,1 9,12 
№8 32678,4 32679,6 2 102,2 60 21,1 9,95 


Чтобы проверить, подчиняются ли данные нормальному закону распределения, построена гистограмма 
частот как один из способов визуального представления распределения данных (рис. 1). Гистограмма имеет 
форму колокола и близка к кривой нормального распределения (рис. 2). 


Король Е.А. и др. Определение сходимости внутрилабораторных результатов измерений запыленности строительной площадки 
и, 


3:5 


3,0 


2,5 


2,0 


Частота 


0,5 


0,0 
7,5 8,9 10,3 И 


Запыленность, мг/м? 


Рис. 1. Гистограмма частот при и = 8 


Рассчитаны и на графике функции плотности нормального распределения приведены среднее значение и 
доверительный интервал проведенных измерений (рис. 2). 


0,35 1,0 
п=8 
0,9 
0,30 в 
0,8 г 
в 0,25 в 0,7 
о [= 
Е Е 
я 
. 0,20 80 
я Е 
№ о 0,5 
5 ы 
© 0,15 Е 
Е Е 0,4 
=] = 
= 5 
Е Е 
0,10 я 0,3 
0,2 
0,05 
0,1 
0,00 0,0 
4 5 6 7 8 9 0. п 12 13 141 15 7 9 14 
Запыленность, мг/м3 Запыленность, мг/м? 
а) 6) 


Рис. 2. Графики нормального распределения при и = 8: 
а — функция плотности нормального распределения; 6 — интегральная функция распределения 


За результат принимается среднее значение с указанием доверительного интервала: 
Киз=9 33 + 0,96 мг/м*. 
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Для определения минимального количества измерений найдены коэффициент вариации, средняя ошибка 
выборки и предельная ошибка выборки (рис. 3). Максимальная статистическая погрешность составила 0,26 мг/мз. 


озоолоовооноовонвовововововосоон000$" —_ ЗОВОВОООООВОВООО000 


® В = 4.15 мг/мз 
10,0 Среднее значение ® 


Кв=0.15 6©=0.49 мг/мМ° и=-1.16 мг/м? 


Запыленность, мг/м? 
© 
> 


0 1 2 3 4 5 6 7 8 9 
Количество измерений, шт. 


Рис. 3. Размах результатов измерений 


Минимальное количество измерений определяется по формуле [11]: 


2.2 

1 
Пит = = , (12) 
Пит =Ь20=2. 


Согласно полученному значению минимального количества испытаний проведены два параллельных 
испытания по определению концентрации пыли на рабочем месте каменщика (таблица 2). 


Таблица 2 
Результаты двух параллельных испытаний 
< Е з г = 
= р = «| & ЕЕ | 98 | вх — 
5 & 5 ЗЕ 8 вч БЕ 5 Е 2 < ах 
5 = И в: [аи аа а=. вЕ. | 8905 ВЕ Зи 
|9 55 я >, оо 2 в я ды = шо в о 
ы |8: 1 З41 28| БОВЕ СНЕ] НЕ 9 а 
ы 5. - Е = ея 8ЕЗ| Ее ВЕ На 
о В: Е | #4 2 ЗЕЯ на ЕЕ ЕВ 
|3 2 | <чк & | ЕВ ЗН | Е ИВ: 
> а [2] баб оо =) я > 
© Е |7) |7) 
п=2 
№1 | 34563,8 | 34565,0 2 102,2 60 22,0 9,98 
9,57 0,59 8,50 8,68 
№2 | 30654,3 | 30655,4 2 102,2 60 22,0 9,15 


Доверительный интервал результатов двух определений концентрации пыли при вероятности 0,95 
составил 0,81 мг/м? (рис. 4). 
0,80 
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0,50 
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0,30 
0,20 
0,10 
0,00 


Плотность вероятности 


Запыленность, мг/м3 


Рис. 4. Функция плотности нормального распределения при и =2 


Король Е.А. и др. Определение сходимости внутрилабораторных результатов измерений запыленности строительной площадки 


Согласно условиям сходимости по ГОСТР ИСО 5725-1-2002 «Точность (правильность и прецизионность) 
методов и результатов измерений. Часть 1. Основные положения и определения» и ГОСТР 51672-2000 
«Метрологическое обеспечение испытаний продукции для целей подтверждения соответствия. Основные 
положения» результаты параллельных измерений можно признать удовлетворительными, так как выполняются 
условия (7) и (10) (рис. 5): 


 =8.68 % < поры = 15 %, 52 =8,50 % < бнзрм = 25%. 

30,0 30,0 
25-4 25,0 
20,0 20,0 


Сходимость результатов, % 
Относительная погрешность, % 


2 3 4 5 6 7 8 
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Сходимость 

— - - Предельный показатель сходимости 

Относительная погрешность 

— - - Максимальная граница относительной погрешности 


Рис. 5. Зависимость сходимости результатов измерений и относительной погрешности 
от количества испытаний 


По результатам обработки принимается: 
Кир =9,57+0,81 мг/м?. 


Определение концентрации пыли весовым методом имеет достаточно высокую составляющую 
погрешностей используемых средств измерений и проводимых операций. Максимально допустимая разница 
результатов двух измерений составляет 0,84 мг/м?. Разница прямых измерений массы двух образцов должна 
быть не более 0,1 мг. При разнице масс в 0,2 мг показатель сходимости результатов двух параллельных 
измерений становится более 18 % и не удовлетворяет условиям прецизионности. При проведении испытаний на 
определение запыленности рассмотренной методикой возможно ограничиться двумя параллельными 
измерениями, при этом погрешность измерений будет равняться инструментальной погрешности. 

Строительная пыль относится к  малоопасным (ТУ класс), разовая предельно допустимая 
концентрация (ПДК) составляет 6 мг/м?, суточная — 10 мг/м3. По результатам испытаний выявлено 
превышение разовой ПДК в 1,6 раза. Для снижения концентрации пыли на рабочем месте каменщика 
рекомендуются комплексные меры по уменьшению пылеобразования на строительной площадке, 
т.к. основными источниками пыли являются подъездные дороги, места выгрузки и погрузки строительных 
материалов, а также технологические процессы, связанные с обработкой и резкой строительных 
материалов [15]. Необходимо проводить мероприятия, снижающие концентрацию пыли на строительной 
площадке и за ее пределами [16]. 

Гигиенический норматив ГН 2.1.6.3492-17 «Предельно допустимые концентрации загрязняющих веществ в 
атмосферном воздухе городских и сельских поселений» устанавливает разовую ПДК строительной пыли в 
городском воздухе не более 0,5 мг/мз. 

Обсуждение и заключение. Определение концентрации пыли весовым методом имеет достаточно высокую 
составляющую погрешностей используемых средств измерений и проводимых операций. Выявлено, что 
максимально допустимая разница результатов двух проведенных испытаний составила 0,84 мг/м3. Разница 
прямых измерений массы двух образцов должна быть не более 0,1 мг. 

Обработка результатов испытаний показала возможность проведения двух измерений по определению 
запыленности рабочего места каменщика на строительной площадке. Сходимость результатов испытаний в 
рамках подтверждения компетентности лаборатории является важным показателем качества проведенной 
оценки и свидетельствует о достоверности полученных данных, поэтому в рамках внутрилабораторных 
сличительных испытаний рекомендуется использовать не менее шести измерений [17]. Чем выше сходимость 
результатов, тем меньше вероятность ошибок и неточностей в оценке условий труда, что, в свою очередь, 
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обеспечивает более объективную информацию о состоянии условий труда на рабочих местах и позволяет 
принять обоснованные меры по их улучшению [18]. 

Работникам исследованного объекта рекомендовано использовать средства индивидуальной защиты органов 
дыхания, защищающие от высокодисперсной пыли. Для снижения концентрации строительной пыли на 
строительной площадке, а также за её пределами следует проводить орошение строительной площадки водой 
для осаждения высокодисперсной пыли, поливать подъездные дороги, промывать колеса автотранспортных 
средств как на въезде на строительную площадку, так и на выезде [19]. 
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Аннотация 

Введение. Вредные компоненты рудной пыли, образующейся при разгрузке продуктов подготовки железорудного 
концентрата (ППЖК) на стадии смешивания, наносят ущерб как работникам, так и оборудованию. Для решения 
этой проблемы применяется орошение жидкостным аэрозолем с использованием форсунок больших диаметров 
(>20 мкм). Однако данный метод неэффективен в улавливании мелкодисперсных частиц пыли, поэтому 
повышение эффективности метода осаждения орошением пыли ПИЖК становится актуальной задачей. Целью 
данного исследования является изучение воздействия технологии «Сухой туман», генерирующей капли жидкости 
размером до 20 мкм, на этапе разгрузки ППЖК горно-металлургического предприятия при осаждении взвешенной 
мелкодисперсной пыли. Основной задачей данного исследования являлась оценка эффективности и возможных 
преимуществ применения технологии «Сухой туман» для орошения пыли с последующим осаждением, поскольку 
к пыли ПИЖК описанная выше технология ранее не применялась. 

Материалы и методы. Эксперимент по осаждению пыли ПИЖК проводился в специально созданном 
лабораторном стенде. Посредством физического моделирования были получены параметры процесса 
осаждения. Далее полученные результаты подвергались анализу с точки зрения получения зависимости 
осаждения пыли с течением времени с учетом влияния технологии «Сухой туман». Для физического 
моделирования была разработана программа эксперимента. Согласно данной программе, пыль равномерно 
загружалась внутрь лабораторного стенда (сверху), распределялась в воздушном потоке по всему объему 
стенда крыльчаткой, а прибор, расположенный в нижней части, фиксировал изменение концентрации во 
времени. Далее были проведены эксперименты по осаждению пыли с применением жидкостного орошения. 
Совместно с загрузкой пыли в объем лабораторного стенда посредством форсунок, генерирующих капли 
размером 10 и 15 мкм, подавалась жидкость — отфильтрованная вода. Эффективность технологии «Сухой 
туман» при осаждении пыли ПИЖК определялась визуально, и далее — на основании сопоставления графиков. 
Изучалась динамика изменения усредненных концентраций пыли от времени как при осаждении без орошения, 
так и с применением технологии «Сухой туман». В процессе эксперимента фиксировались характеристики 
микроклимата внутри лабораторного стенда (влажность, температура и скорость движения воздуха) и 
параметры двух форсунок — их рабочее давление и время распыления подаваемой жидкости. 

Результаты исследования. Сравнение результатов эксперимента показало уменьшение времени осаждения на 
40 % и 75 % при использовании форсунок на 10 мкм и 15 мкм соответственно. 

Обсуждение и заключение. По результатам эксперимента подтверждена эффективность технологии «Сухой 
туман» для осаждения пыли ПИЖК при разгрузке на стадии смешивания. Полученные базисные результаты 
позволят в дальнейшем оценить эффективность осаждения пыли с применением дополнительно режима 
пульсирующей вентиляции. В таком сочетании ожидается повышение эффективности осаждения еще 
на 20-25 % относительно результатов, представленных в данной статье. Полученные результаты дают 
возможность обосновать рациональные параметры и применить на производстве вышеописанный способ для 
повышения эффективности осаждения пыли. Помимо этого, они создают основу для разработки методики 
ускорения осаждения пыли ППЖК с применением метода пульсирующей вентиляции. 
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Введение. В целях снижения концентрации атмосферной пыли на участках производства в 
горно-металлургическом комплексе успешно применяются разнообразные инновационные методы и передовые 
технологии обеспыливания. Используют фильтрационные системы, циклоны и вентиляционные установки, 
специально разработанные для эффективного улавливания и удаления мельчайших частиц пыли [1-3]. Помимо 
этого, инженеры применяют методы жидкостного орошения, аэрофлотации и другие передовые технологии, 
которые совместно обеспечивают комплексный подход к решению проблемы запыленности. Эффективность 
предпринимаемых мер по обеспыливанию тесно связана с уникальными особенностями производственных 
процессов, характеристиками применяемого оборудования, конструктивными особенностями помещений и 
другими важнейшими факторами. Систематическая забота о техническом состоянии оборудования не просто 
является частью технологического процесса, а выступает важным элементом стратегии по снижению 
воздействия производственной деятельности на окружающую среду, способствует оптимизации 
производственной эффективности. Комплексное внедрение вышеперечисленных мероприятий — совмещение 
различных методов и технологий — позволяют достичь максимальной эффективности в уменьшении уровня 
пылевой нагрузки в пределах рабочей зоны. 

На этапе разгрузки и подготовки железорудного концентрата, где формируется значительное количество 
мельчайшей дисперсной пыли, продолжительное время используется система жидкостного орошения [4—6]. 
Для создания мельчайшего жидкостного аэрозоля применяются разнообразные виды форсунок, образующих 
капли диаметром от 30 до 150 мкм. Помимо этого, используются форсунки и с более крупным диаметром. 
Время оседания { мельчайшей дисперсной пыли (4 = 1-10 мкм) при использовании таких форсунок может быть 
достаточно продолжительным и достигать 25 400 с (около 7 часов) [6]. На других этапах производства для 
эффективного пылеосаждения применяется передовая технология — «Сухой туман». В этой технологии 
используются форсунки меньшего диаметра, распыляющие жидкость с дисперсностью капель в диапазоне 
от |1 до 20 мкм [7]. В данной технологии каждая капля жидкости служит эффективным инструментом для 
захвата и обеспечения оседания самых мельчайших частиц, создавая уникальное сочетание технологии и 
инженерного искусства в борьбе с проблемой атмосферной запыленности. 

Исследований, связанных с применением вышеописанной технологии для осаждения мелкодисперсной 
пыли при подготовке продуктов железорудного концентрата, ранее не проводилось. Поэтому целью данной 
работы был анализ воздействия технологии «Сухой туман», при которой генерируются капли жидкости 
размером до 20 мкм, на этапе разгрузки ПИПЖК горно-металлургического предприятия при осаждении 
взвешенной мелкодисперсной пыли. Задачей стала оценка эффективности и возможных преимуществ 
применения технологии «Сухой туман» для орошения пыли с последующим осаждением. 

Материалы и методы. Для оценки эффективности инновационной технологии «Сухой туман» провели серию 
экспериментов пылеосаждения мельчайших частиц продуктов железорудного концентрата с использованием 
форсунок, имеющих диаметры 10 и 15 мкм. Для выполнения эксперимента был разработан лабораторный 
стенд (рис. 1). Данный стенд использовался для анализа осаждения угольной пыли в [8]. Он представляет собой 
кубическую емкость объемом 1 м3, изготовленную из органического стекла в алюминиевой раме. 

В роли средства для мониторинга концентрации частиц в атмосфере использовался высокоточный измеритель 
массовой концентрации аэрозольных частиц — АЭРОКОН-П. Данный измерительный прибор был разработан с 
учетом требований к определению массовой концентрации пыли с разнообразным происхождением и химическим 
составом. Измеритель массовой концентрации обладает уникальной способностью регистрировать дисперсность 
исследуемых частиц, включая те, которые имеют диаметр до 10 мкм. 
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Рис. 1. Лабораторный стенд для осаждения пыли со схемой размещения оборудования и точками измерений: 
1 — генератор воздушного потока; 2 — высокоточный измеритель массовой концентрации аэрозольных 
частиц — «АЭРОКОН-П»; 3 — генератор тумана; 4 — метеометр; 5 — датчик влажности; 6 — датчик температуры; 
7 — форсунка; 8 — засыпная воронка; 9 — воздушный пульсатор 


Порядок проведения измерений обозначены на рис.1 буквами А, Б, В (расположение датчика 
Аэрокон-П) (2). Значения 1, П и Ш (находящиеся рядом с буквами) соответствуют точкам проведения 
измерений. Согласно вышеописанной схеме, эксперимент проводился по всей площади дна лабораторного 
бокса, а именно в 9-и его точках (точки измерения Г, П, Ш следует понимать, как располагаемые в перспективе, 
т.е. друг за другом). 

Эксперименты по осаждению пыли производились при контроле параметров микроклимата. Данные 
параметры на момент проведения эксперимента по осаждению контролировались при помощи метеометра 
ТЕЗ$ТО 435 (4). При проведении экспериментов по осаждению пыли метеометр использовался для контроля 
параметров пылевоздушной среды в лабораторном боксе — температуры и влажности воздуха. Это позволило 
установить начальные условия эксперимента и провести его в контролируемых условиях. Описанный выше 
измерительный прибор был задействован для мониторинга и поддержания влажности воздуха внутри стенда в 
пределах 25—30 % и температуры — в диапазоне от 22 °С до 25 °С. Оценка скорости движения воздушного 
потока, исходящего из генератора, проводилась с использованием цифрового крыльчатого анемометра, что 
обеспечивало стабильную скорость движения воздуха на уровне 4 м/с. Для моделирования системы орошения в 
лабораторном стенде использовалась установка Е218, предназначенная для туманообразования посредством 
распыления жидкости из различных форсунок (рис. 1). Рабочее давление данной установки составляет 5,4 МПа, 
а максимальное — 12,41 МПа. Для эксперимента в данной установке использовались форсунки диаметром 
15 мкм и 10 мкм (рис. 1). 

Размеры частиц пыли, используемой при проведении эксперимента по осаждению, согласно 
гранулометрического анализа, находились в диапазоне от 1 мкм до 40 мкм (рис. 2). Поскольку прибор для 
регистрации концентрации частиц в воздухе измеряет частицы диаметром до 10 мкм, требовалось определить 
процентное содержание этих частиц в пробе пыли. Согласно данному анализу, требуемый размер частиц пыли, 
используемый для проведения эксперимента по осаждению, составляет примерно 10% от 1 грамма навески, 
используемой при анализе пыли [9]. Для корректного определения прибором изменения концентрации навеска 
исследуемой пыли была увеличена до 5 грамм с целью повышения концентрации мелкодисперсной пыли с 
диаметром от 1 мкм до10 мкм. 

Эксперимент по осаждению пыли был реализован в лабораторном стенде, описанном ранее, с соблюдением 
всех параметров микроклимата, регулируемых соответствующими измерительными устройствами. Заданное 
количество пыли массой т=5г вводилось в объем лабораторного бокса в течение 3-5 секунд с 
использованием засыпной воронки (рис. 1). Во время этого процесса генератор воздушного потока был 
активирован, обеспечивая подачу воздуха со скоростью У = 4 м/с. Измерительное устройство АЭРОКОН-П, 
используемое в роли измерителя массовой концентрации аэрозольных частиц, периодически регистрировало 
данные о концентрации пыли с интервалом в 5 секунд, автоматически передавая полученные значения на 
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монитор компьютера [10]. Этап эксперимента по осаждению пыли завершался в момент достижения значения 
концентрации пыли, зарегистрированного прибором, равным п = 0,00 мг/м?. Достижение данного значения 
считалось индикатором завершения эксперимента, при котором весь пылевой материал в лабораторном боксе 


считался осевшим. 
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Рис. 2. Гранулометрический анализ пыли агломерационного производства 
горно-металлургического комплекса: 
С — процентное распределение количества частиц в пробе по размерам, %; 
4С — интегральное распределение процентного количества частиц в пробе по размерам, %; 
4 — размер частиц в пробе 


Используя аналогичные параметры температуры, влажности и скорости движения воздуха при проведении 
эксперимента по пылеосаждению, исследования осаждения пылеводяного аэрозоля проводились при 
поочередном использовании 2-х форсунок (4 = 10 и4 = 15 мкм) и контроле подаваемого давления (р = 5,4 МПа) 
жидкости. Помимо этого, осуществлялся контроль времени распыления {=2 мин (по истечению 2 мин 
влажность воздуха в лабораторном стенде становится максимальной и составляет 98,5 %). 

Результаты исследования. В описанных условиях для обеспечения достоверности экспериментальных 
данных была выполнена серия из 10 экспериментов по осаждению пыли. Результаты данных экспериментов 
были обработаны, сведены в графики и проанализированы. Графики изменения значений концентрации пыли с 
течением времени при проведении 10 экспериментов от времени осаждения приведены на рис. 3. 

В соответствии с графиком, представленном на рис. 3, среднее время осаждения пыли при отсутствии 
орошения и использовании пульсирующей вентиляции составляет в среднем 1 828 секунды (30,5 минуты). В 
течение этого времени пыль в условиях производства находится в воздухе рабочей зоны и наносит вред не 
только оборудованию, но и персоналу предприятия. 

Результаты 10 экспериментов по осаждению пылеводяного аэрозоля с использованием 2-х форсунок также 
были представлены в виде графиков динамики изменения средних значений концентрации от времени 
осаждения аэрозолей (рис. 4, 5). 
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Рис. 3. Графики изменения концентрации пыли от времени по результатам 10 экспериментов 
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Согласно графикам, представленным на рис. 4, среднее значение времени осаждения пылеводяного аэрозоля 
при использовании форсунок диаметром 10 мкм, составляет 1115 с (18,5 мин). Эффективность данного способа, 
относительно времени самоосаждения пыли, составляет порядка 40%. Однако, согласно рис.5, при 
использовании форсунки диаметром 15 мкм, среднее время осаждения пылеводяного аэрозоля составляет 475 с 
(= 8 мин.). Эффективность данного способа составляет 74 %. 
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Рис. 4. Графики изменения концентрации пылеводяного аэрозоля от времени (диметр форсунки 10 мкм) 
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Рис. 5. Графики изменения концентрации пылеводяного аэрозоля от времени (диметр форсунки 15 мкм) 


На рис. 6 представлены пересчитанные зависимости концентрации аэрозолей, полученные в результате 
экспериментов по орошению с использованием форсунок диаметром 10 мкм и 15 мкм. В ходе экспериментов 
были обнаружены разные концентрации пыли, которые были зарегистрированы прибором. Для удобства 
представления эти концентрации были пересчитаны и выражены в виде процентного соотношения, зависящего 
от времени осаждения. 

При сравнении определенных ранее значений эффективности осаждения пыли с применением форсунок 
на 10 и 15 мкм, последние дают значительный эффект. С применением форсунки на 15 мкм значение 
эффективности осаждения пыли примерно в 2 раза выше, чем с форсункой на 10 мкм. При использовании 
форсунки на 15 мкм частицы пыли проявляют низкую склонность к слипанию между собой и менее 
эффективно формируют капли жидкости больших размеров, что снижает эффект «конденсата» в процессе 
орошения. Мелкодисперсные капли лучше заполняют пылевоздушную среду и могут использоваться в тех 
местах технологического производства, где капли больших размеров могут нанести ущерб. 
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Рис. 6. График изменения концентраций пыли (%) от времени в условиях пылеосаждения орошением при применении 
форсунок с диметрами 15 мкм и 10 мкм 


Обсуждение и заключение. В результате проведенных экспериментов по осаждению пыли в лабораторном 
стенде было установлено время осаждения пыли при двух различных сценариях: при естественном осаждении 
пыли и при осаждении пылеводяного аэрозоля (путем воздействия жидкостного аэрозоля с использованием 
форсунок диаметром 10 и 15 мкм). 

В дальнейшем для разработки методики по осаждению пыли необходимо использовать следующие 
параметры, как рациональные: 

— влажность воздуха внутри лабораторного — ф = 25-30 %; 

— температура воздуха в лабораторном стенде — Т= 22—25 °С; 

— скорость движения воздуха, создаваемая генератором воздушного потока, должна быть на уровне /=4 м/с. 

Также следует использовать следующие технические параметры системы орошения: 

— для рассеивания жидкости рекомендуется использовать форсунки с диаметрами 10 мкм и 15 мкм; 

— рабочее давление в системе подачи жидкости — р = 5,4 МПа; 

— время распыления аэрозоля жидкости — 2 мин. 

Полученные результаты показывают разные значения эффективности осаждения пыли с применением 
форсунок, генерирующих разные диаметры капель, применяемых в процессе орошения. Описанные в статье 
форсунки демонстрируют значительное улучшение процесса пылеосаждения орошением при использовании 
форсунок на 10 и 15 мкм, нежели при использовании в процессе орошения форсунок, генерирующих капли 
размером больше 15 мкм [7]. Кроме этого, уменьшение количества влаги в воздушной среде производственных 
цехов будет оказывать меньшее негативное воздействие на оборудование. 

Анализ полученных данных подтвердил более высокую эффективность процесса осаждения пыли при 
использовании технологии «Сухой туман» во время разгрузки ПИЖК на этапе смешивания. Полученные 
результаты предоставляют основу для дальнейшей оценки эффективности процесса осаждения пыли уже при 
использовании метода пульсирующей вентиляции. Ориентируясь на работы других авторов, использующих 
данную технологию [7], прогнозируется дополнительное повышение значений эффективности осаждения на 
20-25 %, относительно применения технологии «Сухой туман». Полученные в данной работе результаты будут 
использованы при разработке в дальнейшем методики осаждения пыли продуктов железорудного концентрата 
с использованием метода пульсирующей вентиляции. 

Полученные значения показателей эффективности технологии «Сухой туман» дают основание 
рекомендовать его для внедрения на производстве для осаждения пыли на участках горно-металлургических 
производств с высокой пылевой нагрузкой. 
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Аннотация 

Введение. Повышенная интенсивность шума и вибрации на железнодорожном транспорте может стать не только 
актуальной проблемой в области охраны труда, но и серьезной экологической проблемой в ближайшем будущем 
из-за необходимости постоянного увеличения веса грузовых поездов, повышения скорости движения 
пассажирских поездов и возрастания пропускной способности железнодорожных линий. В последнее время 
немало сделано для минимизации негативного воздействия шума и вибрации на железнодорожном транспорте. 
Однако проблема превышения показателей виброшумовых факторов остается актуальной. Уменьшение вредного 
воздействия шума и вибрации способствует улучшению условий труда работников поездных бригад и 
повышению комфортных условий для пассажиров. Как правило, основное внимание при анализе воздействия 
виброакустических факторов уделяется внешнему их воздействию на селитебную территорию и 
меньше — воздействию на подвижной состав. К тому же такие исследования не являются комплексными, так как 
не дают полного представления о воздействии шума и вибрации непосредственно на состав в движении. Цель 
данного исследования — анализ воздействия виброакустических факторов на работников поездных бригад и 
пассажиров для разработки метода их комплексной защиты, основанного на применении виброгасящих 
материалов. 

Материалы и методы. В ходе работе были изучены нормативные документы, проведен комплексный анализ 
актуальной информации по данной теме, использованы методики по расчету вибрации и виброускорения. 
Значения уровней звукового давления были получены с помощью шумомера ЗРМ-101. Объектом исследования 
был выбран вагон типа «плацкарт» во время движения. В рамках исследования было измерено звуковое давление 
на выбранных участках дороги и определен класс условий труда. 

Результаты исследований. Результат анализа воздействия виброакустических факторов на работников 
поездных бригад свидетельствует в необходимости усиления комплексных мер их защиты. Полученный расчет 
силы вибрации стал основой предложенного авторами метода минимизации вредного воздействия 
виброакустических факторов, в основе которого использование в конструкции вагона вибро- и 
звукоизолирующих материалов, подходящих по необходимым акустическим параметрам. В частности, 
предложено с целью шумоизоляции на пол вагона уложить плотную резину, а на потолок — материал из 
неорганического волокна. В работе также дана экономическая оценка эффективности мероприятий по 
минимизации вредного виброшумового воздействия. 

Обсуждение и заключение. Предложенные методы минимизации шумового и вибрационного воздействия могут 
способствовать снижению уровня звукового давления и показателей вибрации до нормативных значений, что 
сократит негативное воздействие виброакустических факторов на работников поездной бригады и пассажиров. 
Рассмотренный комплекс по минимизации вибрации и шума включает в себя совокупность методов, в которых 
нашли применение бюджетные материалы, созданные на основе вторично переработанного сырья. 


Ключевые слова: звуковое давление, вибрация, шумовое загрязнение, шумоизоляция 
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Введение. Железнодорожный транспорт всегда вызывал и вызывает шум и вибрацию, которые в большинстве 
случаев не могут быть полностью устранены и, скорее всего, будут негативно влиять на окружающую среду 
невзирая на применение различных методов смягчения их воздействия. Проблема уменьшения влияния шума и 
вибрации железнодорожного транспорта на людей и воздействия на инженерные конструкции очень сложна из- 
за большого количества параметров и факторов распределения шумовой и вибрационной энергии. Негативное 
влияние виброакустических факторов растет с каждым годом в связи с увеличением количества и веса грузовых 
вагонов, большого их износа, несвоевременным техническим обслуживанием пассажирских вагонов, а также 
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нехваткой свободного пространства в городских районах, в результате чего новые деловые и жилые объекты 
размещаются в непосредственной близости от железной дороги. Проблема шума и вибрации от 
железнодорожного транспорта особенно остро стоит в крупных городах с большим железнодорожным 
сообщением главным образом из-за малого расстояния от железнодорожных путей до ближайших зданий. Этот 
факт ограничивает возможности для принятия мер по снижению уровня шума и вибрации. На сегодняшний день 
первостепенной задачей является защита от превышения звукового давления по системе «внешний источник 
шума — окружающая среда». Система подразумевает воздействие шума на близко расположенные объекты 
инфраструктуры, жилые зоны, природу. Защита от такого рода воздействий достигается установкой 
акустических экранов. С их помощью возможно значительно сократить звуковое давление. 

Малоизученной до сих пор является проблема комплексной защиты от виброакустического воздействия. 
Стоит также отметить и недостаток вибродемпфирующих устройств, которые были бы способны снижать 
совокупное воздействие шума и вибрации. По мнению авторов, научная новизна данного исследования 
заключается именно в том, что рассмотрен комплексный подход к решению проблемы минимизации воздействия 
виброакустических факторов и предложены простые шумопоглощающие материалы для массовой их установки. 

Вибрация от рельсового транспорта может оказывать негативные воздействия непосредственно или в 
сочетании с фоновым шумом. Значительное превышение нормативов уровня шума доставляет дискомфорт не 
только работникам и пассажирам, но и жителям прилегающих к железнодорожному полотну территорий, часто 
является значимой причиной ухудшения самочувствия и снижения работоспособности. Специфика шумового 
загрязнения на железнодорожном транспорте такова, что при долгосрочном воздействии в совокупности с 
вибрацией оно очень негативно влияет на работников ж/д транспорта [1]. Длительное и систематическое 
воздействие излишнего шумового давления может привести к развитию хронических заболеваний нервной 
системы, опорно-двигательного аппарата, быть причиной потери бдительности, ухудшения самочувствия, 
ощущения дискомфорта, чувства раздражения [2]. Постоянное воздействие шума даже на номинальных 
значениях (особенно ночью) может стать причиной нарушений сна и хронических неврозов [3]. Количество 
людей с нарушениями сна очень велико — большинство сталкиваются с воздействием постоянного шума в 50- 
60 дБ именно ночью, что при отсутствии других дневных источников шума является одним из основных 
факторов бессонницы [4]. 

Из всех вредных экологических факторов вибрация является одним из самых массовых. Под ее воздействием 
находится от 50 до 70 % населения [5]. Анализ вибраций, возникающих в результате движения подвижного 
состава, очень важен, так как железная дорога — основная транспортная система, охватывающая большую часть 
территории нашей страны и стран ближнего зарубежья. С каждым годом растут объемы и количество 
перевозимых по железной дороге грузов [6]. 

Рабочая вибрация рельсов возникает на стыке колеса и рельса (система «колесо — рельс») и распространяется 
не только через систему опор пути на почвенный покров и окружающие здания, но и непосредственно на 
подвижной состав. В некоторых случаях пассажиры могут непосредственно ощущать вибрацию, которую обычно 
называют вибрацией от взаимодействия подвижного состава и рельсового полотна. Колебания грунта, вызванные 
железнодорожным движением, в основном возникают в зоне контакта между колесом и рельсом [7]. Поэтому 
важно, чтобы действие этого механизма возникновения негативного фактора было сведено к минимуму. Основным 
источником возникновения вибрации на железной дороге является вибрация при качении, вызванная неидеальной 
рабочей поверхностью колеса и дефектами на рабочей поверхности рельса. Эти вибрации передаются колесу и 
путевым конструкциям, что приводит к возникновению превышения пороговых значений по основным октавным 
уровням. Как следствие, вибрация, передаваясь через колесные пары, тележку, вызывает вибрацию кузова вагона, 
что особо ощутимо на высокой скорости поезда. Вибрации высокой амплитуды могут привести к критическому 
состоянию узлов вагона, к урону конструкции рельсового подвижного состава. С целью сохранения оборудования 
вибрации такого рода следует сокращать не только в системе «колесо — рельс», но и непосредственно принимая 
меры по защите основных частей вагона, это, в свою очередь, минимизирует появление вторичной вибрации от 
корпуса вагона. Чтобы свести к минимуму механическое воздействие колеса и рельса, поверхность рельса должна 
быть как можно более гладкой. Вибрации, создаваемые конструкцией вагона и тележки, также должны быть 
сведены к минимуму, чтобы уменьшить количество вибраций, передаваемых на людей в вагоне [8]. Для создания 
системы с желаемыми акустическими и структурно-динамическими свойствами при низком воздействии на 
человека шума и вибрации необходимо учитывать конструкцию пути, колес, тележек, вагонов. Другим важным 
аспектом является обеспечение технического обслуживания этих деталей соответствующим образом 
(своевременное техническое обслуживание и плановый ремонт). 


Соколов Д.А. и др. Обеспечение безопасных условий труда при воздействии на работников поездных бригад ... 


Актуальные методы снижения вибрации, как правило, делятся на две основные категории: категория 
пассивных мероприятий и категория активных мероприятий. В настоящее время чаще всего используются 
пассивные мероприятия. К ним относится применение традиционных виброгасителей, систем амортизации 
ударов разного рода (соударение вагонов) и изоляции основания [5]. Однако современные методы снижения 
вибрации должны быть направлены на пассивно-активное управление вибрацией, которое заключается в 
активном приложении силы, равной противоположным силам, создаваемым внешней вибрацией. Иногда для 
гашения вносимых вибраций предлагается использовать разработанные системы управления и гашения сил, 
которые создают вибрацию, а также конструкции, созданные на основе новых композитных материалов. В 
частности, можно использовать пьезоэлектрические элементы в качестве гашения излишней вибрации'!. Суть 
этого метода заключается в том, что пьезоэлектрический элемент вырабатывает электрическую энергию при 
механическом растяжении и реагирует по форме деформации, если на него подается электрическая энергия. Еще 
одним методом минимизации шума и вибрации является метод применения дисковых тормозов вместо чугунных 
колодок. Альтернативой же более сложным в адаптации, установке и дорогим дисковым тормозам является 
замена чугунных блоков на агломерационные или композитные. Все названные выше методы направлены на 
основные узлы вагона. Авторами был проведен расчет виброизоляции в самом вагоне на пассажирских местах. 

Материалы и методы. В качестве материала исследования был рассмотрен вагон типа «плацкарт». Замеры 
проводились с помощью прибора для измерения звукового давления — шумомера ЗРМ-101. Участки с наиболее 
негативным воздействием виброакустического фактора были выбраны исходя из конструктивных особенностей 
тележки вагона и строения внутренней структуры пространства вагона. Следовательно, план исследования 


включал в себя измерение звукового давления в вагоне и определение класса условий труда. 

Для анализа возможности снижения шума в источнике вагона типа «плацкарт» рассмотрена схема вагона и 
отмечены на ней участки наиболее негативного воздействия виброакустических факторов (рис. 1). На этих 
участках подразумевается наибольшее звуковое давление в вагоне и как следствие — критическое негативное 
воздействие виброакустических факторов при различных условиях (скорость движения, состояние 
железнодорожного полотна и др.). 
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Рис. 1. Значения звукового давления в вагоне по участкам (1, 2, 3) 


Следующим методом исследования минимизации вибрационного воздействия был расчет виброизоляции 
пассажирских мест в вагоне на основе методики?. Были применены виброзащитные вставки из резинового 
материала. Пол вагона у основания кресла вибрирует с частотой /= 50 Гц, с виброскоростью у = 0,4 м/с и массой 
места 22 кг. Примем массу пассажира в 70 кг. Модуль упругости: одо = 0,04 МПа, Е = 2.5 МПа, й=0,1 м, 
коэффициент сопротивления & = 600 Н с/м. 

Определим массу, приходящуюся на вставки, когда пассажир будет сидеть. Масса пассажира, приходящаяся 
на сидение, составляет 70 %, значит: 

т. = 70.0,7 = 49 кг (1) 
тб =49+22 =71кг (2) 


Определим максимальный статический прогиб вставок: 


т 01.4 —16.10-^. (3) 
Е 250 


д 


' Снижение уровня шума в кривых. Железные дороги мира. 2009;6:70-76. ОВГ: В&рз://;4пита.сот/птасез/ра/9т2009-06_70-76.раё (дата 
обращения: 14.11.2023) 

? Пособие по акустической виброизоляции центробежных машин. Москва: Издательство литературы по строительству; 1973. 35 с. 
ОВГ: Врз:/лпегапоги.га/Оаа2/1/4293801/4293801338.р4Ё (дата обращения: 14.11.2023) 
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Частотная характеристика собственных колебаний амортизированного пассажирского места будет: 


в = 2.3.4.65 __ =79с'. 


их 16.10 


о=2л 


Коэффициент передачи вибрации определяется по формуле: 


Найдем относительное демпфирование О и круговую частоту ©: 


Получим: 


= ь _ 600 - 0.05. 
200 -ть 2.79.71 
лы 

90 79 
1+ (2.0,05.4)° 
72 = 0,071. 


[1-4] +(2-0,05-4] 


Затем определим вибрационную скорость на сидении: 


у= 0,4.0,071= 0,028 м/с. 


Частота колебания места составит: 


[0 79 
пе а 


=—= = 12,57 Гц. 
25 2.3,14 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Результаты исследования. Согласно ГОСТу 33787-2019 для частоты 12,57 Гц виброскорость не должна 
превышать 0,0056 м/с3. Полученный в результате расчета показатель не соответствует данному требованию. Для 
его улучшения можно применить виброзащитный материал с другими характеристиками — таким может стать 
композитный лист из плотных резин (рис. 2). Благодаря более высокому коэффициенту сопротивления можно 


добиться минимизации вибрационного воздействия. 


ГОСТ 33787-2019. 


Рис. 2. Вариант расположения виброзащитного материала: 
1 — пассажирское сиденье; 2 — виброзащитная вставка 


Оборудование железнодорожного подвижного состава. Испытания на удар 


Ьрз://4осз.сп4 га/Чоситет 1200170805 (дата обращения: 14.12.2023). 


и 


вибрацию. 


ОВГ: 
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В таблице 1 представлены результаты измерения уровней звукового давления на участках, представленных 
на рис. 1. Сравнив их с нормами уровня звука в вагоне по СП 2.5.3650-20, можно сделать вывод, что наблюдается 
превышение показателей“. 

Таблица 1 
Показатели шумового давления в вагоне по зонам 


Участок, № Уровень звука, дБ 
1 17,5 
2 74,1 
3 78,5 


Проведенные в данной работе расчеты позволили авторам сделать вывод о необходимости применения 
метода минимизации виброакустических факторов, основанного на использовании виброгасящих материалов. 
Его эффективность подтверждается представленным исследованием. Авторами предложен комплекс 
мероприятий, включающий в себя установку в обшивку потолка вагона в наиболее подверженных шуму зонах 
шумопоглощающего материала для минимизации шумового загрязнения. Материал представляет собой 
композит из неорганических волокон толщиной листа от 2 до 4см. Также для снижения шумового и 
вибрационного воздействия предложено установить вставки в пол вагона по аналогии с созданием 
амортизирующей системы. Для расчета эффективности метода был проведен расчет воздействия вибрации на 
пассажирское место с учетом вибрации пола вагона. Материал безопасен при эксплуатации для человека и 
окружающей среды, соответствует критериям пожарной безопасности, которые установлены для пассажирского 
вагона согласно распоряжению ОАО «РЖД» от 5 ноября 2009 г. № 2255р5. 

Обсуждение и заключение. Оценка воздействия вибрации и шума, полученная в данном исследовании, 
показывает, что оно существенно влияет на условия труда работников поездных бригад и пассажиров. Работники, 
которые много времени проводят в дороге, рабочие, занимающиеся ремонтом путей в непосредственной 
близости от железнодорожного полотна, и пассажиры поездов подвергаются крайне вредному шумо- и 
вибрационному воздействию, которое оказывает разрушительное воздействие на их здоровье. Поэтому так важно 
минимизировать это вредное воздействие в источнике его возникновения. Несомненно, улучшение 
существующих и внедрение новых методик минимизации последствий шумового и вибрационного факторов на 
железной дороге положительно повлияет на здоровье работников и будет способствовать созданию 
благоприятной социальной обстановки. Стоит особо подчеркнуть, что мероприятия по минимизации негативного 
воздействия шума и вибрации от железнодорожных составов дают относительно небольшой экономический 
эффект, но имеют при этом большую ценность для обеспечения защиты здоровья как работников, так и 
пассажиров, создания для них более комфортных условий. 

В предложенном авторами комплексе мероприятий по минимизации вибрации и шума применяются 
материалы, созданные из вторично переработанного сырья. Для пола вагона используется плотная резина, для 
потолка — композит из неорганического волокна толщиной 2—4 см. Оснащением вагонов этими материалами 
можно достичь: 

— минимизации шумового давления в вагоне на 7-12 Дб; 

— уменьшения вибрационного воздействия. 


Список литературы 

1. Соколов Д.А. Характеристика шумового загрязнения и вибрации железнодорожного транспорта. 
Минимизация воздействия шума и вибрации на работников железнодорожного транспорта. Студент и наука. 
2022;1(20):65-67. 

2. Капцов В.А., Лексин А.Г. Физические факторы производственной среды и безопасность движения. Гигиена 
и санитария. 2009;5:87-93. ЧВГ: В@рз://субе|йеттКа.гпи/агис1е/п/йлевезК1е-РаКогу-рго1хуод${уеппоу-згеду-1- 
Бехоразпо$-у1яВешуа/\еууег (дата обращения: 14.12.2023). 

3. Епанчинцева Е.А. Экологические проблемы на железнодорожной дороге. В: Материалы 


АТ Международной студенческой научной конференции «Студенческий научный форум — 2019». Москва; 
2019. ОВГ: В&рз://зс1епсеопит.га/2019/агис]е/2018010394 (дата обращения: 10.10.2023). 


* Санитарно-эпидемиологические требования к отдельным видам транспорта и объектам транспортной инфраструктуры. Санитарно- 
эпидемиологические правила СП 2.5.3650-20. ОВГ: № рз://40с$.спё4 ги/4оситеп/ 566406892? тагкег=6560ТО (дата обращения: 14.11.2023). 

5 О введении Инструкции по обеспечению пожарной безопасности в вагонах пассажирских поездов. Распоряжение ОАО «РЖД» № 2255р от 
05.11.2009. ОВГ:  Бйрз:/Леза!ас в ги/дос/тазроцатвеше-оао-г7В4-01-05112009-п-2255г/?узсН4=фау147011074285296 (дата обращения: 
14.11.2023). 


Техносферная безопасность 


В рв:// Брз-оигпа|.га 


Безопасность техногенных и природных систем. 2024;8(1):50-57. е1$5М 2541—9129 


4. Головина Е.И., Соколов Д.А. Обеспечение безопасности труда при эксплуатации железнодорожного 
транспорта путем модернизации систем вибродемпфирующих накладок и рельсовых пластин. Безопасность 
техногенных и природных систем. 2023;(1):39-—46. В рз://401.оге/10.23947/2541-9129-2023-1-39-46 

5. Карабаева М.У., Саидов С.М., Юсуфханов 3.Ю. Распространения вибрации при прохождении поездов в 
зависимости от расположения железнодорожного полотна. Символ науки. 2017;(5):25-27. ЧВГ: В рз://0з- 
газа. соп/ЗВОКМКИЗМ-2017-05.раЁ?узсНа=феКкр9у\9г840350762 (дата обращения: 10.10.2023). 

6. Харламов А.А., Комогорцев М.Г. Метод визуализации и оценки вибрационных воздействий на верхние 
строения железнодорожного пути. Молодой ученый. 2011;3(26):88—91. ПВГ: Бирз://моасВ.га/агстуе/ 
26/2816/?узсПАрэКиеуа2'721587460 (дата обращения: 10.10.2023). 

7. Борисов Е.К., Алимов С.Г., Назаров А.Ю. Вибрация грунта от железнодорожного трафика. Вологдинские 
чтения. 2005;48:87-92. 

8. Романов С.Н. Биологическое действие вибрации и звука: Парадоксы и проблемы ХХ века. Вартанян И.А. 
(ред.). Санкт-Петербург: Наука; 1991. 158с. ОВГ: №@рез://за917фа.сот/е$/зво\и/Фуи/165-1-готпапоу-- 
БоТозтеНезКое-аеуз\1е.В | (дата обращения: 14.12.2023). 


ВегГегепсе$ 

1. боко]оу РА. КВагаК®емз5ЯКа збатоуо2о тазгуахпешуа 1 у1габи тБ@етподого7Впого ‘бтапзройа. 
Мшшихае1уа уо74е1$улуа звита 1 уга5ий па габобКоу тВ@етподого7Вптого ‘гапзрога. 5 и4еп! 1 паика. 
2022;1(20):65—67. (1 Влз$.). 

2. Карбоу УА, Гекзш АС. РВуз1са| Расюг$ оР Фе шдизила| епугоптеп: ап4 1гаЁЯс заРейу. Нуэепе апа Запианоп. 
2009;5:87-93. ПВГ:  Б@рз://субе|ешиКа.га/агас1е/п/Ятасвезке-РаКфогу-рго1иуодуеппоу-згеду-1-6е7оразпоз- 
4у1тВетуа/\1е\уег (ассеззеа: 14.12.2023). (ш Вл$$.). 

3. ЕрапсЬи5еуа ЕА. Епупоптеша! ргоетз оп Фе гаЙ\угау. ш: Маета/[? АТ Мегйаипаго4по! уи4епсйезКо! 
паисрпо Копрегеп5й "Зиаепсйезки паисйпут югит-2019". Мозсо\; 2019. ОВГ: 
6рз://5слепсеогитт.ги/2019/агс]е/2018010394 (ассеззеа: 10.10.2023). (ш Вл$$.). 

4. Сооуша [, ЗоКоюу РА. Епзие оссираНопа| заРефу ай гаЙ\уау тапзрог ореганоп Бу ирогадте Фе зузетв 
оР уШганоп Чатрше рез ап та рез. ба/еу ор Тесйповешс апа Машга! буяет. 2023;(1):39—46. 
0 рз://401.0г/10.23947/2541-9129-2023-1-39-46 (11 Влзз.). 

5. Кагабаеуа МО, За4оу ЗМ, УязиПапоу 7/Уи. Вазргозиапетуа УШгази ри ргокВо7В4епи рое74оу у 
тау15ип1о5Н о газро]о7Нешуа 7Нее7подого7впого рота. бутро[ ор 5сепсе. 2017;(5):25-27. ЧВГ: Б@ре://03- 
газ1а.соп/ЗВОКМИКИЗМ-2017-05.раЁ?узсНЧ=фрэКр9у\9г840350762 (ассеззеа: 10.10.2023). (п Клз$.). 

6. КВа|атоу АА, Котовогзеу Ма. Мею4 у17лаПхазий 1 обзепк1 у1галоппукВ уо74е1$%уй па уегкЬше 
эгоешуа 7т7В@е7тпо4ого7Впо20 рин. Мора исйепу1. 2011;3(26):88—91. ОВГ: 
Брз://тоТасв.ги/агсв1уе/26/2816/?узсИ4=реКиеуа2721587460 (ассеззеа: 10.10.2023). (п Вл$$.). 

7. Вомзоу ЕК, АНтоу $@, Матагоу АУи. Угавтуа огата оё 7е7подого7Впозо 1гаЯКа. Гоюзатяюе сШетуа. 
2005;48:87—92. (ш Вл$5.). 

8. Котапоу 5М. Воюелсйезкое аезйле уфтайй 1 гуика: Рагадокзу 1 ргоЙету ХХ уека. Уацапуап ТА (е4.). 
ЗашЕРеегбиго: МамКа; 1991. 158 р. ОВГ: В@рз://5ба41яба.сот/ез/5во\ми/Чуи/165-1-готапоу--6б1оортевезКое- 
4еузеуе.В | (ассеззеа: 14.12.2023). (ш Вл$.). 


Поступила в редакцию 20.11.2023 
Поступила после рецензирования 16.12.2023 
Принята к публикации 22.12.2023 


Об авторах: 

Дмитрий Алексеевич Соколов, студент кафедры техносферной и пожарной безопасности Воронежского 
государственного технического университета (394006, РФ, г.Воронеж, ул. 20-летия Октября, 84), 
ЗРГ\-код: 3137-6710, ОВС, атипузокоюу598 @етай.сот 


Елена Ивановна Головина, кандидат технических наук, доцент кафедры техносферной и пожарной 
безопасности, заместитель декана по учебной работе факультета инженерных систем и сооружений 
Воронежского государственного технического университета (394006, РФ, г. Воронех, ул. 20-летия Октября, 84), 
ЭРП\-код: 7333-2526, ОВСТО, и00111(@угазч.угп.га 


Соколов Д.А. и др. Обеспечение безопасных условий труда при воздействии на работников поездных бригад ... 


Заявленный вклад соавторов: 

Д.А. Соколов — проведение эксперимента, анализ результатов исследований, подготовка и оформление 
текста статьи, формирование выводов. 

Е.И. Головина — научное руководство, формирование основной концепции, цели и задачи исследования. 


Конфликт интересов: авторы заявляют, что у них нет конфликта интересов. 
Все авторы прочитали и одобрили окончательный вариант рукописи. 


Весеуеа 20.11.2023 
Веукбеа 16.12.2023 
Ассеред 22.12.2023 


АБош фе Аишфогс: 
Отиму А. ЗоКооу, Зва4епЕ оЁ Фе ТесБпозрВеге ап4 Ее За#еёу Оерагитепь Уогопе7В Зе ТесЬт1са! УтуегзИу 
(84, 20 1енуа ОКуабгуа 5+., Уогопе?В, 394006, ВЕ), ЗРТ\-код: 3137-6710, ОВСТО, атинузокоюу598@етай.сот 


Непа Г. Со!оута, Сапа. 5с1. (Еп?.), Аззослайе Ргоеззог оЁ фе ТесБпозрВеге апа Еще Заейу Перагитепь Першу 
Реап Рог Асадепс АНаиз оЁ фе РасиНу ог Епешеегте Зузетз апа Эиасвгез, Уогопе2Н Эа е ТесЬтса! ОшуегзИу 
(84, 20 1епуа ОКуабгуа 5+., Уогопе7В, 394006, КЕ), ЗРПМ-соде: 7333-2526, ОВСТО, и00111(@узази.угп.га 


Сатей соттфшюотуйр: 

РА ЗоКо]оу: сопдисйпе ап ехрегитепь, апау7тс гезеагсЬ гезиЙ$, ргерагте ап4 плате е {ехё оЁ Фе агисе, 
Роге сопс[$101$. 

ЕГ @ооуша: асадепс а4у1$ те, гии айоп оЁ Фе Баз1с сопсерь, 20а! ап4 обуесйуе$ оЁ Фе заду. 


Соп/ПсЕ оГищеге5{ киетет!: ше аийогу 4о по! йауе апу соп1сЕ офпиегези. 


АП аифогу рауе геа4 ап4 арргоуе4 йе дпа тапизсийл. 


Техносферная безопасность 


6 рз:// Брз-оигпа].га 


Безопасность техногенных и природных систем. 2024;8(1):58—72. е1$5М 2541—9129 


ТЕХНОСФЕРНАЯ БЕЗОПАСНОСТЬ 
ТЕСНМОЗРНЕВЕ ЗАКЕТУ 


® не 
ВУ 


УДК 622.8 Научная статья 
6рз://401.0ге/10.23947/2541-9129-2024-8-1-58-72 
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Аннотация 

Введение. Сера присутствует в органической и минеральной части всех видов твердого топлива, причем от ее 
содержания зависит проявление опасных свойств шахтопластов. В связи с этим не теряют актуальность научные 
изыскания, связанные с определением количества серы в ископаемом топливе. Содержание общей серы и ее 
разновидностей изучалось в большинстве случаев с позиций эффективности процессов переработки топлива, качества 
и экологической безопасности угольной продукции. Меньше внимания уделяли проявлению опасных свойств 
шахтопластов, прогнозированию и разработке профилактических мероприятий по предотвращению эндогенных 
пожаров в угольных шахтах. Цель представленной работы — создание инженерного метода расчетного определения 
элементного содержания серы в минеральной массе углей для прогноза опасных свойств шахтопластов. 
Материалы и методы. По результатам корреляционного анализа установили, каким образом содержание серы 
в горючей части топлива зависит от общей серы. Рассмотрели данные по углям практически всех шахтопластов 
Донецкого и Львовско-Волынского бассейнов. Представили соотношения зольности пластовых и обогащенных 
проб, а также сведения о содержании общей серы и серы в горючей части. Это позволило расчетным путем 
оценить содержание серы в минеральных примесях. Исходные данные для разработки метода взяли из 
справочно-нормативных документов. 

Результаты исследования. Рассмотрены семь совокупностей, которые включали от 149 до 1827 шахтопластов. 
По каждой обобщили справочные и расчетные данные: 

— содержание серы — общей (5/7) и в горючей части топлива (5,); 

— условные точки для корректировки показателей во избежание получения отрицательных усредненных значений; 

— эмпирические уравнения, которые описывают расчетные и скорректированные зависимости по соотношению 
5 и 5;; 

— корреляцию этих показателей; 

— диапазоны изменения выхода золы. 

Все эти сведения представлены для обогащенных и не обогащенных углей, для выхода золы более и менее 10 %. 
В одном случае отмечены удовлетворительные результаты определения расчетного содержания серы в условно 
негорючей массе. Это совокупность из 149 шахтопластов. Для 13 из них разность показателей (5/7 — 5.) больше 
или равна нулю. Для 136 — меньше нуля. Для 81 содержание серы в условно негорючей части топлива 5$„ меньше 
нуля. При этом отмечается достаточный уровень корреляции (0,79) и ограниченный диапазон изменения выхода 
золы (2,410). Соответствующее этой совокупности эмпирическое уравнение рекомендуется для прогноза 
опасных свойств шахтопластов с выходом золы менее 10 %. 

Обсуждение и заключение. Полученные результаты можно использовать для совершенствования нормативных 
документов по безопасному ведению горных работ с учетом содержания серы, как в органической массе, так и в 
минеральных примесях ископаемого угля. 
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Введение. Сера — элемент соединений, входящих в органическую и минеральную часть любого топлива. В 
минеральной массе углей сера встречается главным образом в виде сульфатов (в основном кальция и железа) и 
дисульфидов железа Ре5> (пирита и марказита), в органической массе — в виде органических соединений. 

Различают несколько основных форм серы в углях. Это части общей серы: 

— сульфатная 550, — входит в состав неорганической массы угля в виде сульфатов металлов; 


— пиритная 5, — входит в состав неорганической массы угля в виде дисульфидов металлов (пирита и марказита); 
— органическая 5, — входит в состав органических соединений. 
Общую серу и ее разновидности, в отличие от серы золы, условно пересчитывают на элементную серу [1]: 

5: =550, +5р +90. (1) 


Содержание серы в золе принято рассчитывать в виде триоксида серы (503). Сера — нежелательная и вредная 
часть топлива. При сжигании угля она выделяется в виде 5О», загрязняет окружающую среду, вызывает коррозию 
металлических поверхностей, уменьшает теплоту сгорания топлива. При коксовании сера переходит в кокс, 
ухудшая его свойства и качество металла. По указанным причинам содержание общей серы и ее разновидностей 
изучалось в большинстве случаев с позиций эффективности процессов переработки топлива, качества и 
экологической безопасности угольной продукции. В этом направлении получены хорошие результаты, что 
подтверждается разработкой и успешным использованием современной промышленной классификации углей по 
генетическим и технологическим параметрам. В значительной степени она базируется на данных справочников, 
в которых по итогам экспериментов обобщалась информация о качестве каменных углей и антрацитов [2], а 
также об их обогатимости [3]. 

Менее значимы успехи в прогнозировании и разработке профилактических мероприятий по предотвращению 
эндогенных пожаров при ведении работ в угольных шахтах. Так, в течение двадцати лет (с 1994 по 2015 год), в 
шахтах Украины произошло более двухсот эндогенных пожаров [4]. Согласно открытым данным, аналогичные 
аварии происходят во многих угледобывающих странах. Таким образом, не решена проблема, актуальная с 
ХХ века. Вначале она была связана с самовозгоранием угля на складах и при транспортировке, затем стали 
происходить аварии в угольных шахтах, обусловленные эндогенными пожарами [5]. 

Содержание серы в угле не только определяет его потребительские качества, но и существенно влияет на 
проявление опасных свойств шахтопластов при ведении горных работ [6]. В частности, при увеличении 
содержания в угле общей серы на 1—6 % вероятность эндогенного пожара возрастает более чем в три раза [7]. 
Были попытки разобраться, как по минеральным примесям можно выявить опасные свойства угольных 
шахтопластов [8]. Однако и сейчас почти не известно, каким образом исследуемые риски зависят от вида и 
формы серы в ископаемом угле. В результате в нормативных документах"? вероятность эндогенных пожаров в 
шахтопластах определялась только содержанием общей серы?. 

Исследования [5] показали, что самовозгорание в значительной мере обусловлено природой не только угля, но 
и сопровождающих пород“. Крайне важно установить составную часть топлива, от которой зависит вероятность 
самовозгорания. Без этого лабораторные методы исследований не дают приемлемых ответов на поставленные 
вопросы. Предполагалось наряду с лабораторными исследованиями угольного вещества учитывать наличие в 
ископаемом топливе углистой породы из минеральных пропластков, а также из пород кровли и почвы. Для этого 
необходимо дополнительно рассматривать выход золы и содержание всех видов серы [5]. 

В руководстве по предупреждению и тушению эндогенных пожаров°, как и в других нормативных 
документах по безопасному ведению горных работ, не принимается во внимание наличие минеральных 
примесей. Они не учитываются при борьбе с пылью в угольных шахтах, прогнозах газовыделения° и внезапных 
выбросов? и в целом в нормативных документах, регламентирующих горные работы. 


' Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000 от 18.12.2000. Донецк: 
Научно-исследовательский институт горноспасательного дела; 2000. 216 с. 

2? Каталог углей СССР, склонных к самовозгоранию. Москва: Недра; 1982. 416 с. 

3 Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000. Донецк: НИИГД; 2000. 216 с. 
* ГОСТ25543-2013. Угли бурые, каменные и антрациты. Классификация по генетическим и технологическим параметрам. ОВГ: 
В ре://е.55оушЕ/а2/1/4293772/4293772639.раЁ (дата обращения: 11.01.2023). 

5 Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000 от 18.12.2000. Донецк: 
Научно-исследовательский институт горноспасательного дела; 2000. 216 с. 

° Руководство по проектированию вентиляции угольных шахт. Киев: Основа; 1994. 311 с. 

7 Правила ведення арничих роб?т на пластах, схильних до газодиналчних явищ. Стандарт Мнвуглепрому Укра!ни № СОУ 10.1.00174088.011- 
2005. Кив: Мвуглепром Украйни; 2005. 224 с. 
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Содержание минеральных примесей, судя по выходу золы из пластовых проб угля [2] и [3], может 
значительно превышать 40 %. Это, несомненно, сказывается на достоверности прогноза опасных свойств 
шахтопластов при ведении горных работ и указывает на актуальность совершенствования нормативной базы. 
Для этого необходимо разработать инженерный метод определения серы в минеральной массе углей. 

Цель представленной работы — доказать принципиальную и практически единственную возможность 
расчетного определения элементного содержания серы в минеральной массе углей для прогноза опасных свойств 
шахтопластов при ведении горных работ. Отметим, что при реализации данного подхода применяются некоторые 
стандартные показатели, зафиксированные в справочно-нормативной документации по потребительским 
качествам топлива для каждого шахтопласта. Речь идет об экспериментальных данных, касающихся содержания 
серы — общей и в горючей части топлива. Благодаря этим сведениям можно установить соотношения между 
содержанием серы в органической массе и минеральных примесях. По выходу золы можно определить 
соотношения между органической массой и минеральными примесями для обогащенных проб угля. 
Соотношение зольности пластовых и обогащенных проб дает возможность оценить содержание серы в 
удаляемых продуктах обогащения. 

Материалы и методы. Методика проведения исследований основана на тесной корреляционной связи 
("= 0,93) содержания серы в горючей части топлива (5,) и обшей серы (5/2) для углей практически всех 
шахтопластов Донецкого и Львовско-Волынского бассейнов (рис. 1). 


5% 


5% 
0 1 2 в. 4 5 6 ре: 


Рис. 1. Зависимость содержания серы в горючей части топлива (5.) от общей серы (5/)) согласно статистической обработке 
экспериментальных данных: 1 — усредняющая прямая; 
2 — биссектриса координатной сетки [2,3] 


Статистическую обработку зависимости 5; от 5/А выполнили для 1827 шахтопластов, по которым известно 
содержание серы — общей и в горючей части топлива. 

В справочном нормативном документе [3] указаны общие сведения о потребительских качествах углей для 
734 шахтопластов, а в справочнике [2] — для 1379. 

Близость усредняющей прямой | и биссектрисы 2 координатной сетки (рис. 1) показывает, что в общем 
случае содержание серы в горючей части топлива (5:) может быть выше и ниже значения общей серы (5/). 
Разница между 5 и 5, (5 -5,) в первом приближении может быть критерием соотношения элементного 
содержания серы в горючей массе и в условно негорючей части топлива. Показатель (5-5) 
для 1315 шахтопластов был отрицательным, что указывает на преимущественное содержание серы в горючей 
части массы. Для 429 шахтопластов у этого показателя были положительные значения, что свидетельствует о 
преимущественном нахождении серы в условно негорючей части топлива. Нулевая разница (5 - 5,) 
для 83 шахтопластов означала равенство содержания серы в горючей и условно негорючей частях углей для всей 
совокупности шахтопластов, рассматриваемых в статье. 

Имеется информация для 2113 шахтопластов. По некоторым из них совпадают данные в разных документах, 
и это аргумент в пользу достоверности сведений. Для четырех шахтопластов не приведены значения выхода золы 
из пластовых проб (таблица 1). По этой причине можно анализировать экспериментальные данные о выходе золы 
для 2109 шахтопластов. Из них обогащение проводилось для 1622 шахтопластов, а для 487 не проводилось. 
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Таблица 1 


Сведения о шахтопластах, для которых не приведены данные о выходе золы из пластовых проб 


Страница 


Индекс 


Зольность, % 


Содержание серы, % 


Источник Шахта пластовой обогащенной - са | Вгорючей 
источника пласта общей, 5; 
пробы, Аи пробы, Ао° части, 5. 
№ 144 Г — 6,3 6,3 1,9 
[3] 74 В Аз = 6,1 2,8 3,0 
10-бис тз — 7,0 3,2 3,5 
[4] 106 Им. Войкова Ь = —- — _ 


Результаты исследования. Показатели из [2] и [3] позволили рассчитать для 1827 шахтопластов содержание 


серы в условно негорючей части топлива 5: 


5 = 


Н 


100-54 — (100-45) -5, 


где Ао’ — выход золы из обогащенных проб, %. 
Расчетные значения 5, связаны с показателем (5/_5,) тесной прямо пропорциональной корреляционной 


зависимостью (т = 0,91, рис. 2 а). 


5 


н? 


% 


А“ 


[9] 


6) 


‚% 


5 


ЭГ, % 


55, % 
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Рис. 2. Зависимость от критерия (51 — 5:) расчетных значений содержания серы (5н) по уравнению (2) в условно негорючей 


части топлива: а — для всей совокупности из 1827 шахтопластов; 6 — для 1352 шахтопластов, угли которых 
предварительно обогащались; в — для 475 шахтопластов с необогащенными углями; 1 — усредняющие прямые; 
1' — скорректированные положения усредняющих прямых; А' — условное нахождение точки для корректировки 
расположения прямой 1' и исключения получения отрицательных («диких») усредненных значений 5н; Хх — расчетные 
значения 5н по уравнению (2) [3, 4] 


Параллельно с этим получено значительное количество отрицательных расчетных значений 5, для 
802 шахтопластов. Тесная корреляция показателей и наличие «диких» (отрицательных) результатов 
свидетельствуют о возможном влиянии систематических погрешностей при определении разницы между 
содержанием серы общей (5/') и в горючей части топлива (5,). Наличие отрицательных расчетных значений 5, 
может указывать на то, что систематические погрешности определения критерия (5/- 5.) обусловлены как 
занижением экспериментальных значений общей серы, так и завышением в горючей части топлива. 

Для раскрытия природы появления «диких» расчетных значений содержания серы в условно негорючей части 
топлива по уравнению (2) рассмотрели методики определения показателей, входящих в эту зависимость. 

Известны несколько подходов к определению общей серы в углях — от классического гравиметрического 
метода Эшка до современных видов инструментального анализа [1]. Метод Эшка, предложенный в 1874 году, 
все еще широко используется в мировой практике. Он регламентируется ГОСТом. Навеску угля сжигают путем 
спекания со смесью оксида магния и безводного углекислого натрия (соотношение 2:1 по массе) в окислительной 
атмосфере при температуре (800 + 250) °С. Органическая масса топлива сгорает, и сера превращается в сульфаты 


натрия и магния. Их переводят в раствор. Сульфат-ион количественно осаждают в солянокислой среде 
хлористым барием в виде сернокислого бария Ваз Ол. По массе Ва5 О. рассчитывают общее содержание серы. 
Недостатки этого подхода: большие затраты времени (не менее 6-7 часов) и возможность потери части 
органической серы при спекании со смесью [1]. В большинстве случаев метод применяется при возникновении 
разногласий как арбитражный. Общую серу в стандартных образцах топлива, предназначенных для градуировки 
анализаторов, определяют только методом Эшка. 

Для практических целей рекомендуется регламентированный ГОСТом? ускоренный стандартный метод 
определения общей серы. Навеска топлива сжигается в потоке кислорода или воздуха в трубчатой печи при 
1350 °С. Образуются зольный остаток и газообразные продукты сжигания, в состав которых входят оксиды серы 
(в основном 502) и хлора. Анализ выполняется за 40-60 минут. Этот метод использовали при определении 


$ Правила ведення арничих робт на пластах, схильних до газодинащчних явищ. Стандарт Минвуглепрому Украни СОУ 10.1.00174088.011- 
2005. Кив: Манвуглепром Украйни; 2005. 224 с. 

ГОСТ 2059-95 (ИСО 351-96). Топливо твердое минеральное. Метод определения общей серы сжиганием при высокой температуре. Москва: 
Стандартинформ; 2008. 14 с. ОВГ: 6 рз://Шез.зноушЁ.та/Оайа2/1/4294832/4294832840.р4{ (дата обращения: 11.11.2023). 
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содержания общей серы для установления потребительских качеств углей, экспериментальные значения 
приведены в справочно-нормативных документах [2,3]. Ускоренный метод определения общей серы 
рекомендуется для контроля за содержанием серы в малозольных топливах. Его нежелательно использовать при 
анализе высокозольных топлив со значительным содержанием соединений кальция. Из-за разнообразия состава 
минеральной массы и соединений серы нельзя быть полностью уверенным в том, что часть серы не останется в 
золе в виде СабОл. Как следствие, результаты получаются заниженными [1]. 

Для большинства шахтопластов после обогащения угля выход золы составляет менее 10 %, но этот показатель 
бывает и значительно выше (рис. 3). 
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Рис. 3. Зависимость выхода золы из обогащенных проб угля (Ао°) от выхода золы из пластовых проб (Акс) согласно 
экспериментальным данным [2, 3] 


По этой причине даже после обогащения значительной части проб нет оснований считать все угли малозольными. 
Для 487 шахтопластов, согласно данным [2, 3], обогащение не проводилось и выход золы во многих случаях превышал 
10%, что также не позволяет относить угли этой совокупности шахтопластов к малозольным. 

Соединения кальция есть в минеральных примесях углей всех шахтопластов. Об этом свидетельствует 
типизация состава зол (таблица 2), согласно которой содержание СаО в общем случае может достигать 40 %. 


Таблица 2 
Типизация состава зол углей по данным [2] и [3] 
. Пределы колебания состава, % 
Типы зол углей з 

АБОз УХ @ Ее2Оз Сао 

Кремнистые 8—30 40—70 до 20 до 20 

Глиноземные 30—45 40—55 до 20 до 20 

Железистые 10-20 30-55 >20 до 20 

Известковые 5—20 15—40 5—20 20—40 
Одна из причин получения «диких» расчетных значений 5, по уравнению (2) — занижение 5 при 


определении ускоренным методом. Такой вывод обосновывается: 

— анализом выхода золы из обогащенных и пластовых проб; 

— возможным содержанием соединений кальция во всех типах зол. 

Второй фактор получения отрицательных расчетных значений 5, — завышение доли горючей части топлива 
с использованием в качестве множителя (100—А.°) в уравнении (2). Содержание общей серы 5/ относится как к 
органической части, так и к минеральной массе (ММ) топлива. По этой причине более логичным было бы 
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использовать в уравнении (2) вместо выхода золы Ао’ значений ММ, которые для рассматриваемых шахтопластов 
экспериментально не определялись. Согласно известным эмпирическим зависимостям, значения ММ всегда 
превышают выход золы на сухое состояние. Учет такого превышения требует дополнительных (недостающих) 
сведений о гидратной влаге, выделяющемся диоксиде углерода при полном разложении карбонатов, пиритной и 
сульфатной сере. Для приблизительной оценки содержания минеральной массы иногда используют соотношение 
ММ = 1,1 4о* [1]. 

Анализ показал причины систематического занижения произведения 100-5/ и завышения разности 
(100-/5°) 5: при определении расчетных значений 5, по уравнению (2). 

Наличие отрицательных значений 5, для совокупности из 802 шахтопластов (рис. 2 а) дает возможность 


оценить среднее занижение критерия в. — к ‚› При котором получены «дикие» результаты: 


Ш 


25-5) 
т. 


0,42, % (3) 


Здесь п = 802 — количество шахтопластов, для которых получены отрицательные значения 5; 2(5* —5. ) а 
1 


сумма отрицательных значений 5/ — 5, для совокупности из и = 802 шахтопластов. 


Среднему отрицательному значению р —5, =- 0,42% соответствует среднее отрицательное значение тм : 
п 
2.5. 1682,18 
Я в = —=-2,10,%, (4) 
п 802 
п 
где р 5; — сумма отрицательных расчетных значений 5, для совокупности из п = 802 шахтопластов. 


1 


Чтобы исключить получение отрицательных усредненных значений нь необходимо положение прямой 1 


(рис. 2 а) скорректировать с учетом координат к — 5, =-0,42 % и 5 =—-2,10 %. Скорректированная прямая 1' 


должна проходить не ниже точки А' с координатами о В, =-0,42% и 5, =+2,10%. Прямо 


пропорциональная зависимость 5„ от (5-—5,) (г=0,91) указывает на то, что истинное положение 
скорректированной прямой 1' будет близко к ее параллельности с прямой 1. Использование в инженерных 
расчетах значений 5,„', определенных по уравнению 1' (рис. 2 а), гарантирует, как минимум, отсутствие их 
отрицательных усредненных значений. 

Эмпирическое уравнение 1" получено на основании экспериментальных данных для 1827 шахтопластов. Из 
них для проб углей 1352 шахтопластов проводилось предварительное обогащение. Выход золы из пластовых 
проб углей остальных 475 шахтопластов был несколько ниже, поэтому для них обогащение не 
предусматривалось. 

Предварительное обогащение существенно повлияло на выход золы по сравнению с необогащенными 
пробами (рис. 3). Следовало выяснить, как точность расчетного определения серы по уравнению (2) в условно 
негорючей части топлива зависит от обогащения. Для этого рассмотрели отдельные характерные совокупности 
шахтопластов. К одной совокупности из 1352 вариант отнесли шахтопласты, угли которых подвергались 
предварительному обогащению (рис. 2 6). До обогащения выход золы из пластовых проб составлял 11,0-48,2 %, 
а после — 3,7—16,7 % (рис. 3). 

Во второй совокупности из 475 вариант рассмотрели шахтопласты, угли которых предварительно не 
обогащались (рис. 2 в), так как условно их относили к малозольным. Выход же золы из необогащенных углей 
этой совокупности шахтопластов колебался от 2,4 % до 28,1 % (рис. 3). 

Выход золы для обеих рассматриваемых совокупностей шахтопластов во многих случаях значительно 
превышал 10 %. Диапазоны выхода составляли 3,7—16,7 % и 2,4-28,| % соответственно. Это указывает на 
отсутствие в ряде случаев для обеих совокупностей обоснованного применения альтернативного ускоренного 
стандартного метода определения общей серы 57 [9]. В результате по большинству рассматриваемых 
шахтопластов расчетные значения 5/ оказались заниженными. Вследствие этого получены и «дикие», 
рассчитанные по уравнению (2), значения 5, соответственно для 558 и 244 шахтопластов обеих (2-й и 3-Й) 
совокупностей (таблица 3). 
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Средние отрицательные значения 58 —5 для шахтопластов 2-й и 3-й совокупностей были рассчитаны по 


уравнению (3). Они равнялись соответственно —0,42 % и —0,41 %. Это указывает на практически одинаковые 
средние погрешности занижения расчетного содержания общей серы альтернативным ускоренным методом для 
обеих совокупностей при расчете 5, по уравнению (2). По указанным обстоятельствам ускоренный метод 
определения обшей серы и возможные занижения ее истинных значений не относятся к различиям между 2-й и 
3-й совокупностями шахтопластов. 

Логика уравнения (2) показывает, что разные расчетные значения 5 2-й и 3-й совокупностей можно получить 
только с учетом разного выхода золы из обогащенных и пластовых проб углей. Для удобства анализа 
потенциального влияния выхода золы из обогащенных и необогащенных проб углей уравнение (2) представили 
в виде слагаемых. 

Для обогащенных проб (совокупность 2): 


1 1 
5+5. (5) 


Для необогащенных пластовых проб (совокупность 3): 


п |] 1 | п 
5 ОО а (6) 
А® А т т 


п п 


Сравнивая слагаемые уравнений (5) и (6), определили направленность их влияния на расчетные значения $5,. 
а а 
Первое слагаемое (100/ А ) .5, уравнения (5) заведомо болыше аналогичного слагаемого (10с/ А )-5. 


уравнения (6). Это определили по бесспорному превышению выхода золы из пластовых проб по сравнению с 
выходом из обогащенных (Ан > Ао°). Неравенство позволяет утверждать, что значения общей серы, определенные 
ускоренным альтернативным методом, будут в большей степени занижаться для пластовых проб. По этой 
причине значения 5/ уравнения (5), как минимум, не меньше, чем в (6). 


Вторые слагаемые со знаком минус уравнений (5) и (6) ((100/4$).5° и (100/4=)-5"') характеризуют 


увеличение доли содержания серы в горючей части топлива после обогащения углей по сравнению с 
содержанием в пластовых пробах до обогащения. Под воздействием процессов обогащения выход золы 
сокращается в несколько раз (4о° < Аи‘, рис. 3). В результате увеличивается доля горючей части топлива 


обогащенных углей (100/42) по сравнению с ее долей в пластовых пробах (100/4° ). Это, а также 


устойчивость соединений серы обусловливают рост относительного (процентного) содержания серы Ку: 
горючей части топлива после обогащения. Эти значения всегда превышают процентное содержание серы в 
пластовых необогащенных пробах, т. е. соблюдается неравенство 9: > м. Отрицательные слагаемые 


уравнений (5) и (6) ((100/4$).5° и (100/4=)-5") при любом их соотношении составляют лишь некоторую 


часть от содержания серы в горючей массе. Вследствие этого суммарные значения второго и третьего 
слагаемых уравнений (5) и (6) будут положительными. 


Рассмотрим неравенство слагаемых в (5) и (6): 5° > 5" и (100/45) 5% и (100/4)-5%. Это говорит о том, 


что обогащение ведет к процентному увеличению содержания серы в негорючей части топлива (5н) и ее 
тесной, почти функциональной корреляционной зависимости (т = 0,94) от критерия 5/5, (рис. 2 6). Сильная 
корреляция бн и 5! — 5, для обогащенных углей шахтопластов 2-й совокупности обусловлена значительным 
принудительным сокращением содержания условно негорючих (минеральных) примесей. Об этом говорит 
сокращение диапазона выхода золы из пластовых необогащенных проб с 11,0-—48,2 % до 3,7-—16,7 % после 
обогащения. Для рассматриваемой совокупности шахтопластов в большинстве случаев (около 90 %) выход 
золы составлял менее 10 %, что указывает на незначительную долю содержания минеральных примесей в 
углях после обогащения. Обогащением в значительной степени объясняется приближение массовой доли 
горючей части топлива к органической массе. Корреляционная зависимость | (рис. 2 6) относится к 
незначительной доле (менее 10 %) сохранившихся после обогащения минеральных (условно негорючих) 
примесей. Относительно высокий коэффициент регрессии уравнения 1 (13,15) указывает на диспропорцию 
показателей: 

— сокращения выхода золы после обогащения углей; 

— содержания серы в сохранившейся части условно негорючих (минеральных) примесей [10]. 
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Как видно из уравнения (5), высокая корреляционная зависимость 5» от показателя 51 — 5, для шахтопластов 
2-й совокупности обеспечивалась относительным постоянством выхода золы Ао° для большинства шахтопластов 
в пределах 3,7—10,0 %. Это фактически исключает существенное влияние непредсказуемых колебаний доли 
условно негорючих (минеральных) примесей в более широких диапазонах. Значения 5„° для 2-й совокупности 
шахтопластов, рассчитанные по уравнению (5), определяются экспериментальными данными 5/ и 55. Это 
подтвердили результаты статистической обработки (рис. 2 6). Влияние колебаний выхода золы Ао’ после 
обогащения углей свелось до минимума, и расчетные значения 5„° в большей степени зависели от двух 
показателей — 5Ё и 5,5. 

Такой устойчивой связи между 5, и (5/-5,) не наблюдается при отсутствии процессов 
обогащения (рис. 2 в). Это выразилось в более низком коэффициенте корреляции (7 = 0,56) и снижении 
коэффициентов регрессии эмпирических уравнений 1 с 13,15 (рис. 2 6) до 3,50 (рис. 2 в). 

При расчете значений 5," по уравнению (6) для 3-Й совокупности шахтопластов диапазон 
непредсказуемых колебаний значений выхода золы Ап расширился до 2,4 %-28,1 %. Это значительно 
больше, чем при стабильном выходе золы из обогащенных проб углей 2-й совокупности. Отметим усиление 
влияния дополнительного фактора, связанного с существенными колебаниями выхода золы из пластовых 
проб в более широких интервалах. Это привело к ослаблению корреляционной зависимости | до значения 
т = 0,56 (рис. 2 в). 

Подавляющее большинство углей 3-й совокупности шахтопластов по критерию 54—-5,<0 должно 
относиться к малозольным. Из 475 шахтопластов всей рассматриваемой совокупности для 454 (т. е. для 95,6 %) 
значения критерия были отрицательными. По отрицательным значениям критерия 5 — 5, основное процентное 
содержание серы необходимо было условно отнести к горючей массе углей этих шахтопластов. Определение 
5 согласно уравнению (6) для 244 шахтопластов дало «дикие» значения. Для остальных 231 шахтопласта 
расчетные значения 5н были положительными (210 вариант) или равнялись нулю (21 варианта). Исследование 
выявило примерно равное количество околонулевых расчетных значений 5„. Это подтверждает наличие 
систематических погрешностей экспериментально определенных параметров, входящих в (6). 
Доказательством служит также примерное равенство по абсолютной величине средних отрицательных 
(-0,71 %) и средних положительных (+0,80 %) расчетных значений 5,. Это дает возможность обоснованно 
вносить количественные корректировки в расчеты, чтобы исключить получение «диких» результатов. 

Сравнение расчетных значений содержания серы в негорючей части топлива для 2-й и 3-й совокупностей 
демонстрирует значительное отличие исходных параметров показателей, которые использовались при 
определении 5, из," по уравнениям (5) и (6) соответственно. Основная разница обусловлена 
предварительным обогащением углей 2-й совокупности. Обогащение непропорционально изменило 
соотношение выхода золы и содержания серы в горючей части топлива. Исходные соотношения 
формировались в природных условиях и связаны с индивидуальной геологической историей образования 
каждого шахтопласта, а именно: 

— последовательным накоплением осадков; 

— погружением и преобразованием исходного вещества на предшествующих метаморфизму стадиях 
торфообразования и диагенеза. 

В углях 3-й совокупности шахтопластов сохранились исходные сложившиеся природные соотношения между 
выходом золы из пластовых проб, содержанием общей серы и некоторой ее части в горючей массе топлива. 

Обогащение углей 2-й совокупности шахтопластов и сокращение выхода золы (в большинстве случаев 
ниже 10%) предопределили искусственную близость всей совокупности по содержанию минеральных 
(условно негорючих) примесей. В результате фактического устранения влияния выхода золы из 
обогащенных проб для большинства шахтопластов по уравнению (5) получены расчетные значения 5н", 
которые в значительной степени непосредственно связаны с критерием (5 - 5,). Это и предопределило их 
высокую корреляционную взаимозависимость (рис. 2 6). 

Главное отличие углей 2-й и 3-й совокупностей шахтопластов — значения выхода золы из обогащенных 
и пластовых проб углей. Для детального представления о влиянии выхода золы из обогащенных и 
необогащенных проб углей на расчетные значения 5„ дополнительно рассмотрели четыре совокупности 
шахтопластов с углями разной степени глубины обогащения (рис. 4). 


Пронская Н.В. и др. Определение содержания серы в минеральной массе для прогноза опасных свойств угольных шахтопластов 


5% | 5% 


545, % 


154-5., % 


-1,0 -0,5 0,0 0,5 1,0 2,0 3,0 -2,0 -1,5 -1,0 -0,5 0,0 0,5 1,0 


в) 2) 
Рис. 4. Зависимость расчетных значений содержания серы (5н) в условно негорючей части топлива от критерия (5 — 5, [2, 3] 


По выходу золы из обогащенных проб шахтопласты разделили на две группы: менее 10 % и более 10% 
(рис. 3). Рассмотрели соответствующие группам совокупности — 4-ю (рис. 4 а) и 5-ю (рис. 4 6). 

Большинство шахтопластов (1181 % или 87,4 %) по выходу золы менее 10 % из обогащенных проб отнесли к 
4-й совокупности. В 5-й совокупности шахтопластов с выходом золы более 10% рассмотрели расчетные 
значения 5 для 171 шахтопласта, что составило 12,6 % от общего количества шахтопластов (1352), угли которых 
обогащались. 

Для обеих совокупностей шахтопластов после обогащения получены тесные прямо пропорциональные 
корреляционные зависимости расчетных значений 5, от критерия 5/ — 5,. Коэффициенты корреляции — 0,95 
и 0,92 соответственно. Коэффициенты регрессии этих зависимостей (рис. 4а, 6) равнялись 13,84 и 7,50 
соответственно. Почти двукратная разница между эмпирическими коэффициентами регрессии свидетельствует 
о существенных различиях между углями 4-Й и 5-й совокупностей. Ее можно объяснить: 

— количественным и качественным составом минеральных примесей в пластовых пробах; 

— разной эффективностью обогащения. 

Влияние указанных факторов требует отдельного детального исследования [11]. 

Аналогичные совокупности шахтопластов сформировали по выходу золы из пластовых (необогащенных) 
проб менее и более 10 %. Для 6-й совокупности из 149 шахтопластов установили зависимость 5н от критерия 
5—5; при выходе золы из необогащенных (пластовых) проб менее 10 % (рис. 4 в). В 7-й совокупности из 326 
шахтопластов рассмотрели аналогичные зависимости расчетных значений 5н от критерия (5 — 5,) при выходе 
золы более 10 % (рис. 4 г). 

Полученные прямо пропорциональные эмпирические зависимости для 6-й и 7-й совокупностей шахтопластов 
с разным выходом золы из необогащенных (пластовых) проб значительно различаются. Отличия обусловлены 
как теснотой корреляционных связей, так и коэффициентами регрессии эмпирических уравнений. 
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Выход золы из необогащенных (пластовых) проб менее 10 % указывает на низкую зольность углей 6-й 
совокупности шахтопластов. В этом случае корреляционная связь довольно высока (г = 0,79). Для более зольных 
углей 7-й совокупности шахтопластов корреляция значительно слабее (у = 0,56). 

Более чем в три раза различаются коэффициенты регрессии эмпирических уравнений 6-Й и 7-й совокупности 
шахтопластов. Они соответственно равны 8,27 (рис. 4 в) и 2,53 (рис. 4 г). 

Выявленные существенные отличия между эмпирическими зависимостями 6-й и 7-й совокупностей шахтопластов 
по тесноте корреляционных связей и коэффициентов регрессии указывают на разные особенности условий 
накопления серы в минеральной (условно негорючей) части топлива для углей с разным уровнем выхода золы. 

Низкий уровень выхода золы (менее 10 %) свидетельствует о незначительном содержании минеральных 
примесей. Между расчетным содержанием серы в условно негорючей части и критерием (5/ — 5,) (рис. 4 в) 
фиксируется довольно высокий коэффициент корреляции (г= 0,79). Это значит, что сера формировалась 
параллельно в минеральных примесях и в органической массе на разных стадиях образования угля. Отсутствие 
такой тесной корреляции для высокозольных углей указывает на разные условия формирования содержания серы 
в минеральных примесях и в органической массе 7-й совокупности шахтопластов (рис. 4 г). 

Обсуждение и заключение. Согласно эмпирическим уравнениям (рис. 4) попарно рассмотрели совокупности 
шахтопластов, которые при прочих равных условиях отличались по двум факторам (таблица 3): 

— обогащались или не обогащались угли; 

— выход золы из обогащенных или необогащенных углей был менее или более 10 %. 

Обогащение существенно меняет соотношения между: 

— органической массой и минеральными примесями; 

— горючими и условно негорючими частями топлива. 

По этой причине 1-я, 2-я, 4-я и 5-я совокупности шахтопластов и соответствующие им эмпирические 
уравнения не отражают истинное природное содержание серы в негорючих частях ископаемых углей. Их нельзя 
использовать для прогноза опасных свойств шахтопластов при ведении горных работ. 

Содержание серы в условно негорючей массе необогащенных проб для 3-Й и 7-й совокупностей занижается 
из-за наличия высокого выхода золы и применения альтернативного стандартного метода определения 
содержания общей серы. Область применения соответствующих этим совокупностям эмпирических уравнений 
ограничена шахтопластами с низкозольными углями и незначительным содержанием соединений кальция в 
минеральных примесях. 

Наиболее достоверные результаты определения расчетного содержания серы в условно негорючей массе 
получены для 6-й совокупности шахтопластов. Практическое применение соответствующего эмпирического 
уравнения рекомендуется при прогнозе опасных свойств шахтопластов с выходом золы менее 10 %. 

Можно улучшить предложенный метод определения серы в минеральных примесях и расширить область его 
применения. Для этого необходимо установить соотношения: 

— содержания серы в органической и горючей массе углей; 

— минеральных примесей и выхода золы (с определением ее химического состава). 

Результаты таких исследований можно использовать для совершенствования нормативов по безопасному 
ведению горных работ с учетом содержания серы как в органической массе, так и в минеральных примесях 
ископаемых углей. 
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Аннотация 


Введение. Вопросы обеспечения оптимальных режимов расхода топлива являются предметом постоянного 
исследования в области эксплуатации автотранспортных средств. Работы в области снижения расхода топлива 
в большинстве случаев посвящены исследованию конструкции силовых установок, трансмиссии или 
аэродинамических особенностей кузова автомобилей. Вместе с тем вопросы определения оптимальных законов 
управления неустановившимся движением с точки зрения синтеза алгоритмов управления движением 
автомобиля не рассматривались. Для оптимизации передаточного числа автомобилей при неустановившимся 
движении авторами работы проведены аналитические расчеты, позволяющие упростить синтез алгоритмов 
управления движением, что позволяет снизить расход топлива. Целью работы являлось определение 
оптимального передаточного числа, позволяющего решить задачу синтеза управления транспортным средством 
для снижения расхода топлива. 

Материалы и методы. Рассматриваемый в работе синтез алгоритмов управления движением автомобиля 
основан на применении к инвариантным признакам действительного движения игольчатой вариации 
Л.С. Понтрягина. Использован аналитический метод оценки энергетической эффективности эксплуатационных 
характеристик автомобиля, в основе которого лежит определение оптимального передаточного числа 
трансмиссии автотранспортных средств с учетом наименьшего расхода топлива. Представленный метод 
учитывает величину крутящего момента, передаваемого от коленчатого вала двигателя на элементы трансмиссии, 
который, в зависимости от мощности двигателя, реализуется в виде силы тяги на колесах автомобиля. 
Результаты исследования. Построен закон оптимального изменения передаточного числа трансмиссии при 
разгоне автомобиля за минимальное время. Решена задача определения оптимального передаточного числа 
трансмиссии автомобиля в случае движения автомобиля с постоянной скоростью и постоянной подачей 
топлива, а также при разгоне автомобиля до заданной скорости при постоянной подаче топлива, когда 
соблюдается условие = = соп$. Результатом рассмотренного случая применения оптимального закона 
изменения передаточного числа трансмиссии является минимизация расхода топлива при ограничениях на 
ускорение (силу тяги) и скорость движения автомобиля. 

Обсуждение и заключение. Применение авторами аналитического метода для определения передаточного числа 
трансмиссии автомобиля, а также использование данного метода в практических расчетах для автомобиля с 
заданными характеристиками показывает возможность решения задачи синтеза управления транспортным 
средством с применением математического аппарата. Это подтверждается построенной графической 
зависимостью по результатам проведенных расчётов. Рассмотренные случаи движения позволили определить 
аналитические зависимости оптимального передаточного числа и скорости автомобиля. Исходные данные, 
полученные по аналитическим зависимостям, применимы для автомобилей со смешанным режимом управления. 


Ключевые слова: передаточное число трансмиссии, расход топлива, энергетическая эффективность 
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Введение. При движении автомобильного транспорта, с учетом его перемещения в условиях городских 
застроек, режим управления транспортным средством является смешанным с чередованием разгона и 
торможения, а также с необходимостью различного маневрирования. Подобное управление приводит к 
повышенному расходу топлива. Уменьшение расхода топлива в Таких условиях, как правило, связано с 
необходимостью увеличения времени совершения маневров [1]. Это вызывает необходимость поиска 
оптимальных законов управления неустановившимся движением автомобиля. В работе [2] энергетические 
характеристики перемещения транспортного средства определяются в первую очередь величиной 
передаваемого от коленчатого вала двигателя крутящего момента и реализуемой на движителях силы тяги, 
которые зависят от мощности двигателя и значения передаточных чисел. Необходимость определения 
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оптимальных значений передаточных чисел трансмиссии транспортного средства обусловлена уменьшением 
значений топливных характеристик [3], а также уменьшением величины энергии, которая тратится при 
маневрировании автомобиля в случае неустановившегося движения. Передаточные числа определяются в 
соответствии с заданной характеристикой двигателя и в зависимости от скорости перемещения транспортного 
средства, имеющей переменные значения [4]. 

При рассмотрении применяемых подходов для решения задач, связанных с оптимизацией, следует сказать 
об их многообразии. Большая часть таких подходов использует методы оптимального управления и 
вариационного исчисления. Авторы в [5] использовали принцип максимума Л.С. Понтрягина с определением 
аналитических зависимостей для гибридных автомобилей. Подобный подход рассматривается и в работе [6], на 
основе которого предлагается построение закона оптимального изменения передаточного числа трансмиссии 
при разгоне автомобиля за минимальное время, а также определение топливной экономичности в ездовом 
цикле [7]. Анализ проведенных исследований показал значительный теоретический задел в области 
определения оптимального передаточного числа при различных режимах эксплуатации автомобиля. Вместе с 
тем, следует сказать, что оптимальный закон изменения передаточного числа трансмиссии должен 
определяться с учетом минимизации расхода топлива при ограничениях на ускорение (силу тяги) и скорость 
движения автомобиля. Авторами представленной работы предлагается для определения оптимальных законов 
управления неустановившимся движением (с точки зрения синтеза алгоритмов управления) применить метод 
оптимизации, в основе которого лежит применение к инвариантным признакам действительного движения 
игольчатой вариации Л.С. Понтрягина. Предлагаемый подход авторами ранее проведенных исследований не 
рассматривался. 

Целью работы является определение оптимального передаточного числа, позволяющего решить задачу 
синтеза управления транспортным средством для оптимизации расхода топлива. 

Материалы и методы. При расчете мгновенного расхода топлива автомобиля используем следующую 
зависимость: 

Ст =ЧуКоКмМ., (1) 
где 4х — показатель удельного расхода топлива, соответствующий максимальному значению № = Миах; 
К» — значение коэффициента, определяемого с учетом скоростного режима работы двигателя; Ку — значение 
коэффициента, определяемого с учетом степени загрузки двигателя; № — мощность двигателя в текущий 
момент времени. 

Коэффициенты Кь и Км определяли с учетом эмпирических зависимостей, К» — с учетом скоростного 
режима работы двигателя. Определим его по формуле: 
2 


о’ ь › (2) 
у 9 


где а» = 1,27; Бь = 0,94; с. = 0,67; о. — текущее значение угловой скорости коленчатого вала; «ом — угловая 
скорость вращения коленчатого вала, соответствующая максимальной мощности двигателя № нах. 
Величина @. связана с передаточным числом трансмиссии { и скоростью движения автомобиля Хх 


зависимостью: 
в, =(Г.х)/г.. (3) 
Тогда формула (2) примет вид: 
Ь , с РЕ 
К, =а,-—® (я) (1-4). (4) 
мт, (®мг, 


Коэффициент Кх зависит от степени загрузки двигателя 5. Значение = определяли отношением мощности 
двигателя в текущий момент времени независимо от режима движения транспортного средства. Принимали В 
расчет: текущее значение угловой скорости коленчатого вала @е мощность двигателя № в момент полной 
подачи топлива и такое же значение угловой скорости коленчатого вала ©е: 


=—°. (5) 


Значение текущей мощности двигателя № определяли дифференциальным уравнением движения 
автомобиля по оси х, имеющего вид [5]: 


о. М : 
Вере ур. (6) 
8 тах 

Здесь Ге) —_ значение коэффициента, определяемого [© учетом затрат МОЩНОСТИ С учетом приращения 


кинетической энергии вращающихся масс двигателя, трансмиссии и колес при разгоне автомобиля; 
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\’ — коэффициент полезного действия трансмиссии; т — полная масса автомобиля; 
лу = Г со5а + та — коэффициент сопротивления движению; К — коэффициент обтекаемости, определяющий 
силу сопротивления воздуха. 

Коэффициент 5 в представленной работе определили по следующей зависимости б = 1 +7, в которой 
значение у зависит от конструктивных параметров автомобиля и изменяется в пределах у = 0,001-0,003. Тогда 
из уравнения (6) можно получить выражение для текущего значения мощности двигателя: 

_ тах | 1+& ь 

Пт 8 

Значение мощности двигателя №. в момент максимальной подачи топливной смеси определили с помощью 
известной зависимости (формула Лейдермана) для внешней характеристики двигателя, имеющей вид: 

2 
т: аа Е (8) 
м ® м © 


М 


Е: 


чи (+ |. (7) 


0) 
= М, тах 


М 


ес 


Эмпирические коэффициенты 4, В, С, зависящие от параметров двигателя, могут определяться, например, 
по зависимостям, приведенным в работе [5]. После преобразований с учетом (3) выражение (8) принимает вид 


№, =Е(ай+ЬРх? -ср хо), (9) 


где коэффициенты а, В и с определили по формулам: 


М В С 
Е". ДЕД. р ВИННИ 
7 ВР Ом — (ох 
Эмпирическая зависимость для коэффициента К» имеет вид: 
Ков ее а (10) 
где константы ам = 3,27; м = - 8,22; см = - 9,13; ам = 3,18 для бензиновых двигателей; ам = 1,20; Вх = 0,14; 
см = 1,80; 4м = — 1,46 для дизельных двигателей. 


Степень загрузки двигателя = в зависимости (10) с учетом (5), (7) и (9) определили по формуле: 
.2 

те®уи, | 1+ .. : 

ВИ ну (ИР 

М тахЙт 8 

= 2 (р 
и+Ьгх-сгх 

В качестве целевой функции использовали величину расхода топлива за время 1», в течение которого 


осуществляется движение автомобиля на траектории х(1)<х, с переменной скоростью х= { (Е). Задача 


оптимизации решали в соответствии со следующим алгоритмом: устанавливали закон, в соответствии с 
которым определяется изменение передаточного числа трансмиссии Их) › И соответствующую ему траекторию 


х(И с таким учетом, чтобы целевой функционал принимал минимальное значение: 


к 
Л = | бе > тт, (12) 
0 
передаточное число, текущие скорость и мощность удовлетворяли ограничениям: 
2, ХХ, М, < Ми, (13) 


где и — передаточное число главной передачи. 
Результаты исследования. Результаты проведенных расчетов были использованы для случая 
непрерывного изменения передаточного числа трансмиссии. Тогда, пренебрегая инерцией вращающихся масс 


(при у = 0), оптимальное значение ик, если оно существует, находится из условия: 
д(@тх) (ок, _дКн К 
О ОК К РАО (14) 
6 & —@ Ку 6 


Полученное выражение (15) в качестве результатов исследования применили для двух предельных случаев. 


Первый случай — движение автомобиля с постоянной скоростью Хх = соп5{ и постоянной подачей топлива, 
когда 2 и Км = соп5{. Поэтому второе слагаемое в уравнении (14) равно нулю, что позволяет получить простую 
формулу для определения оптимального передаточного числа трансмиссии: 

—- Бом 7, к 


Ею 15 
ыы 26,х ео 
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При минимизации времени разгона до заданной скорости оптимальное передаточное число трансмиссии 
определили по формуле [5]: 


Вог ЗАС 
М, =- Ик 1+ |+ (16) 
"ЗС: В? 
Второй случай —Щ_ разгон автомобиля до заданной скорости при постоянной подаче топлива, когда 


соблюдается условие = = сопз(. В этом случае значение коэффициента К» = соп5{ и оптимальное передаточное 
число трансмиссии определили по формуле из решения уравнения: 


ОК. х 
х+Ко —- 0. 
0} 0: 
После преобразований с учетом выражения (4) для К», получим: 
Ь Ь х 
бр на, 9 р) (1-0 (17) 
Фи, (вк. ) мг. (фу, ) 0 
0х 


Производную Ра в уравнении (17) определили из условия & = соп5Е с использованием выражения (11). 
й 


При этом будем пренебрегать, как и ранее, инерцией вращающихся масс (у = 0) и силой сопротивления воздуха 
(кг = о) . Тогда связь между скоростью движения автомобиля и передаточным числом трансмиссии установим 
(при & = сой5й) из решения следующего уравнения: 
За о . 2 
СВ? -ЬР-а+ (#) 
54 


= 0, 
где введены обозначения: 


% ы 
= тт р) = чу (и |. 
т м7, 8 


Полученное уравнение относительно скорости имеет решение, если его дискриминант больше или равен нулю. 


3 
4с 


а 
Тогда для передаточного числа трансмиссии должно выполняться условие: 


Р=Ь? +4асй — 


2(1)> 0. 


ее 18 
"О ы(В? /4С+А) — 


Соответствующее этому условию значение скорости определили по формуле: 


При разгоне автомобиля предельное (минимальное) значение передаточного числа трансмиссии Тит 


0х 
определили из условия (18), что позволяет получить достаточно простое выражение для производной Ру и 
[ 


распространить его на всю траекторию движения: 
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Подставляя полученные выражения в уравнение (17), после преобразований получим уравнение для 
оптимального значения передаточного числа трансмиссии: 
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Для примера провели расчет расхода топлива для автомобиля со следующими характеристиками: полная 
масса автомобиля т = 1800 кг; коэффициенты внешней характеристики двигателя А = 0,64; Б = 1,36; С = 1,0; 


ит = 


максимальная мощность двигателя нах = 100 кВт; угловая скорость вращения коленчатого вала, 
соответствующая максимальной мощности двигателя ол = 576 с'!; радиус колеса 7% = 0,34 м; КПД трансмиссии 


Птр = 0,85; передаточное число главной передачи 1, =3,4 и передаточные числа коробки переключения 


передач 15 = 1,0; 4 = 1,5; в=2,2; 2 =3,2; п =4,8; коэффициент вращающихся масс 6 = 1,0; коэффициент 
сопротивления движению \у = 0,12. 

Приведенный пример расчета расхода топлива для автомобиля с заданными характеристиками наглядно 
показывает решение задачи синтеза управления транспортным средством. Графики разгона автомобиля со 
ступенчатой коробкой переключения передач, иллюстрирующие решением для оптимального закона изменения 
передаточного числа, приведены на рис. 1. Для сравнения представлен график для разгона транспортного средства 
со ступенчатой коробкой переключения передач, который был построен по полученной аналитической зависимости. 


На первом этапе разгона до скорости х„„ кривые совпадают. В данном случае х/ин = 6.14 м/с или 22,1 км/ч. 
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Рис. 1. Графики разгона автомобиля со ступенчатой коробкой переключения 
передач (сухая укатанная дорога а = 100; /= 0,025; у = 0,118): 
орг — оптимальное передаточное число трансмиссии; 
«т — передаточное число, определяемое по времени разгона 
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ГДе Флтах — коэффициент мощности двигателя; Г», — начальная скорость машины на 1{-той передаче; 
1 =1,3 — индекс, обозначающий номер передачи; /х0) — скорость машины на /-той передаче. 

Обсуждение и заключение. Применение авторами аналитического метода для определения передаточного 
числа трансмиссии, а также использование данного метода в практических расчетах для автомобиля с 
заданными характеристиками показывает возможность решения задачи синтеза управления транспортным 
средством с применением математического аппарата. 

Поставленная авторами работы цель — определение оптимального передаточного числа, позволяющего 
решить задачу синтеза управления транспортным средством для снижения расхода топлива — была достигнута. 
Результатом проведенных исследований являются аналитические зависимости, полученные для двух 
предельных случаев движения автомобиля. 

Проведенные авторами расчеты для первого случая показали оптимальное значение передаточного числа 
трансмиссии, которое определили в соответствии с формулами (15) и (16). Значение передаточного числа 
трансмиссии обратно пропорционально скорости движения автомобиля. Значение скорости, определяемое для 
оптимального передаточного числа трансмиссии будет меньше, чем передаточное число, определяемое по 
времени разгона. Формирование эмпирической зависимости во втором случае при переменной подаче топлива 
Ки = К® и оптимальное передаточное число трансмиссии были определены в результате численного решения 
уравнения (14). 

Важность проведенного исследования заключается в упрощении решения задачи синтеза управления 
неустановившимся движением транспортного средства, что наглядно показывают полученные графические 
зависимости. Подобный результат возможен при применении объединенного принципа максимума. 
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Содержание и дисперсность ферросплавов в обмазке 
при микродуговом легировании стали 


М.С. Степанова, Ю.М. Домбровский ® 
Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация : 
зберапоутз(@уапаех.га ЕРМ: МРЕИВ 


Аннотация 

Введение. Основным недостатком традиционных процессов диффузионного поверхностного упрочнения 
стальных изделий является большая продолжительность, поэтому проблема интенсификации таких процессов 
является актуальной. Для ее решения предложено применение высокоэнергетического воздействия на материал, 
позволяющего получить упрочненный поверхностный слой из обмазки, в состав которой входят порошки 
ферросплавов, содержащие легирующие элементы. В литературе отсутствуют данные о необходимом 
содержании и дисперсности таких порошков в составе обмазки. Цель исследования — выбор размера частиц 
ферросплавов и их концентрации в обмазке для достижения максимально эффективного упрочнения 
обрабатываемого изделия. 

Материалы и методы. Для экспериментальных исследований использовали цилиндрические образцы из стали 
20 диаметром 12мм и длиной 35 мм, на поверхность которых наносили легирующую обмазку, содержащую 
порошки ферросплавов и электропроводный гель в качестве связующего. После этого образцы погружали 
вертикально на половину длины в металлический контейнер, который далее заполняли угольным порошком с 
размером частиц 0,4-0,6 мм. Затем пропускали электрический ток величиной от 2,5 до 3,0 А в цепи источник 
питания — контейнер — угольный порошок — образец. Продолжительность процесса составляла 2—8 мин. 
Результаты исследования. Выполнена расчетная оценка электрической проводимости угольного порошка, 
рассчитаны теплофизические параметры микродугового нагрева стали: мощность, выделяемая электрическим 
током на поверхности стального изделия, плотность теплового потока, энергия единичного микродугового 
разряда. Получены выражения для расчета размера частиц порошка ферросплавов, а также экспериментальные 
зависимости толщины диффузионного слоя от размера частиц ферросплавов и их содержания в обмазке. 
Обсуждение и заключение. Результаты исследования позволили определить диапазон размеров частиц 
ферросплавов и их содержание в обмазке, позволяющие обеспечить наиболее эффективное легирование 
поверхности обрабатываемых изделий. Полученные данные будут использованы при разработке 
технологических процессов поверхностного упрочнения стальных изделий. 


Ключевые слова: химико-термическая обработка, энергия микродуги, диффузионное насыщение стали 
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Введение. Продолжительность различных процессов химико-термической обработки (ХТО) стали может 
достигать 8—10 часов [1,2]. Поэтому ранее были предложены методы их интенсификации, основанные на 
применении высокоэнергетического воздействия на материал с использованием плазменного нагрева [3], 
электро-химико-термической обработки [4], комбинированных способов нагрева [5], лазерной обработки [6], 
микродугового оксидирования [7], нагрева в электролите [8]. Однако эти методы сложны для практического 
применения и поэтому не получили широкого распространения. Для интенсификации процессов ХТО 
целесообразно применение микродугового поверхностного легирования, которое позволяет достичь 
значительного (в десятки раз) ускорения процесса диффузионного насыщения, и не требует использования 
сложного оборудования. Особенностью данного метода является погружение обрабатываемого изделия с 
нанесенной на его поверхность обмазкой в угольный порошок с последующим пропусканием через них 
электрического тока. Интенсификация процесса насыщения достигается за счет образования ионов диффузанта 
и их последующей ускоренной диффузии в материал под действием электрического и температурного полей. В 
ранее проведенных исследованиях достигнуты положительные результаты применения данного метода для 
диффузионного насыщения хромом, молибденом, ванадием [9], вольфрамом [10]. Для комплексного насыщения 
на поверхность обрабатываемого изделия предварительно наносилась обмазка, в состав которой входил порошок 
ферросплава, содержащего диффузант. Эти исследования позволили определить структуру и фазовый состав 
полученных покрытий, но задача достижения максимальной эффективности упрочнения, для решения которой 
необходимо определение оптимальных размеров частиц ферросплава и его содержания в обмазке, в ранее 
проведенных исследованиях не ставилась. 
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Цель исследования — определение условий максимально эффективного упрочнения стальных изделий 
методом микродугового легирования за счет выбора размеров частиц порошков ферросплавов и их концентрации 
в используемой обмазке. 

Материалы и методы. Для экспериментальных исследований использовали образцы из стали 20, содержащей 
0,2 мас. % С, диаметром 12мм, на поверхность которых наносили обмазку, приготовленную на основе 
электропроводного геля с добавкой порошков ферросплавов. Образцы погружали вертикально на половину длины в 
металлический контейнер диаметром 35мм, который далее засыпали угольным порошком с размером частиц 
0,4—0,6 мм. Далее через образец пропускали электрический ток величиной от 2,5 до 3,0 А в течение 2-8 мин. 

Для приготовления обмазки использовали порошки феррохрома, ферромолибдена, феррованадия, 
ферровольфрама с размером частиц 0,40-0,50 мкм и связующее в объемном соотношении 1:1. 

Образцы подвергали шлифовке и полировке по стандартной методике с последующим травлением реактивом 
Ржешотарского. Микроструктуру диффузионного слоя исследовали с помощью микроскопа Меорво-21. 

Для измерения насыпной плотности угольного порошка использовали методику по ГОСТ 32558—2013. 
Порошок рассеивали по размерам частиц с помощью набора сит по ГОСТ 33 029-2014. 

Результаты исследований. Насыпная плотность угольного порошка, т.е. масса единицы объема свободно 
насыпанного порошка определена экспериментально, полученные значения приведены на рис. 1. 


П, г/смз 
0,90 т 


0,70 = = - х, мм 
0,2 0,4 0,6 0,8 1,0 


Рис. 1. Зависимость насыпной плотности П угольного порошка от размеров Х его частиц 


В диапазоне размеров от 0,4 до 0,6 мм среднее значение насыпной плотности можно принять равным 
0,81 г/см3. Полученное значение позволило определить объемную плотность угольного порошка как отношение 
насыпной плотности к плотности каменного угля, значение которого заимствовано из справочной 
литературы [11] и принято равным 1,6 г/смз, откуда объемная плотность принята равной 0,81/1,6 = 0,51. 
Полученное значение примерно соответствует плотности упаковки простой кубической решетки (0,52). 
Поэтому, согласно [12], можно принять, что общее сопротивление порошковой среды представляет 
собой сумму сопротивлений последовательных слоев, состоящих из параллельных цепочек из 
контактного сопротивления Ак между частицами и сопротивление Ам самих частиц (рис. 2). 


Рис. 2. Электрические контакты угольных частиц в порошковой среде: 


Кк— сопротивление межчастичных контактов, Юм — сопротивление отдельной угольной частицы 


Расчет по данным рис. 2 позволил оценить общее электрическое сопротивление К порошковой среды: 


РЕ у» +Аи 4” 

— 21 Н(+а7) 

где т — количество слоев угольных частиц, / — радиус образца; 4 — диаметр частиц, Н — глубина погружения 
образца, Ак — сопротивление межчастичного контакта; Ам — сопротивление отдельной угольной частицы. Для 
расчетов были приняты следующие значения: 4 = 0,5103 м, Н = 1,5:10- м, г = 6 мм, Ак = 6 кОм, Ам = 0. Значение 
В принято равным 15,5 Ом. Полученное значение позволило рассчитать: мощность электрического тока Р = РВ, 


(0 
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энергию единичного микродугового разряда Обь, поверхностную плотность тока /, плотность 4 теплового потока. 
Величина 0 определена как отношение мощности тока к суммарному количеству микродуг, которые в течение 
1 секунды одновременно воздействуют на нагреваемую поверхность. По данным видеозаписи значение М 


принято равным 22 500. Результаты расчетов приведены в таблице 1. 


Таблица 1 
Значения теплофизических параметров микродугового нагрева 
ТА [.102, А/см? Р.10?, Вт а: 10°, Вт/м? О-103, Дж 
2,50 45 0,97 1,71 4,3 
2.75 49 117 2.07 52 
3,0 53 1,39 2,47 6,2 


Далее рассчитали максимальный размер Чтах частицы ферросплава, ДЛЯ парообразования которой достаточно 


энергии одной микродуги. 
Необходимое количество теплоты О) с учетом удельной теплоты сублимации /, можно записать как: О = Ат, 


где т — масса частицы. 


Для сферической частицы плотностью р: 


т=рИ = р .. 
6 
Поэтому искомое условие имеет вид: 
Арта? /6< 0). 
Максимальный диаметр частицы, необходимый для выполнения этого неравенства: 
Я ых ЕЗ м ь (2) 
пр 
Рассчитанные значения 4 шах представлены в таблице 2. 
Таблица 2 
Результаты расчетов 
Ферросплав Ее+Ст Ее+Мо Ее+у Ее+ 7 
0,45 А/см? 0,045 0,043 0,044 0,038 
Ялах, ММ 0,49 А/см? 0,048 0,046 0,047 0,040 
0,53 А/см? 0,050 0,048 0,050 0,042 


Таким образом, установлено, что максимальный размер частиц ферросплавов, используемых в составе 
обмазки, не должен превышать 38 мкм. 

Экспериментально определена толщина диффузионного слоя в зависимости от объемной доли частиц феррохрома 
в обмазке и ее толщины. Для этого в связующее добавляли порошок феррохрома ФХО10А с размером частиц 
40—50 мкм в количестве, необходимом для получения его объемного содержания в пределах 10—50 %. 

Результаты представлены на рис. 3. Видно, что максимальная толщина покрытия образуется при объемной 
концентрации частиц феррохрома в обмазке 50% и ее толщине 0,5 мм. Таким образом, максимальная 
интенсификация процесса диффузионного насыщения достигается при условии наибольшего содержания частиц 
ферросплавов в обмазке, а воздействию микродуг должны подвергаться частицы ферросплава, прилегающие к 
насыщаемой поверхности. 


Н, мкм 


100; 


50 — 


0,5 1,0 1,5 2.0 


Рис. 3. Зависимость глубины Н упрочненного слоя от толщины В слоя обмазки для различного содержания порошка 
феррохрома (об. %) в обмазке: 1 — 10;2 — 30; 3 — 50 
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Обсуждение и заключение. Выполненная расчетная оценка электрической проводимости угольного 
порошка, используемого для микродугового легирования стальных изделий, позволила установить влияние 
величины электрического тока в цепи на теплофизические параметры процесса: плотность теплового потока на 
нагреваемой поверхности, мощность, выделяемую на нагреваемой поверхности, энергию единичного 
микродугового разряда, возникающего между нагреваемым изделием и окружающей порошковой средой при 
протекании электрического тока. Полученные результаты позволили определить диаметр частиц ферросплавов, 
используемых для микродугового легирования, их объемную долю в обмазке, которые обеспечивают достижение 
наибольшей интенсификации процесса микродугового легирования. Результаты исследований будут 
использованы при разработке технологических процессов поверхностного упрочнения стальных изделий 
методом микродугового поверхностного легирования. 
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Стимуляция бейнитного сценария превращения внешним 
магнитным полем 
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Аннотация 

Введение. Изменение свойств сталей с бейнитной структурой имеет практический смысл, т. к. при бейнитном 
превращении под воздействием магнитного поля возможно улучшение пластичности стали при сохранении или 
повышении ее прочностных показателей. Научные изыскания в этой сфере касались вопросов влияния 
магнитного поля на термодинамику и смену сценария фазового превращения. Однако в открытых источниках 
нет детального описания воздействия магнитного поля на структуру и свойства продуктов промежуточного 
бейнитного превращения. Цель работы — исследование особенности влияния внешнего магнитного поля на 
сценарий и кинетику фазового превращения стали. 

Материалы и методы. Исследование проводилось на образцах из стали 65Г. Их химический состав 
контролировали при помощи оптико-эмиссионного спектрометра МазеПап 08. Термическую обработку 
(резистивный нагрев) проводили в установке для высокотемпературных исследований «ИМАШ 20-75». 
Температура нагрева — около 1000 °С, время выдержки — 10 минут. Образец охлаждали при помощи 
водоохлаждаемых электроконтактов. Внешнее магнитное поле напряженностью 400кА/м и 800 кА/м 
создавалось электромагнитом, интегрированным в вакуумную камеру установки. 

Результаты исследования. Эксперименты подтвердили возможность смены сценария превращения с 
перлитного на бейнитный при воздействии внешним магнитным полем до 1 МА/м. Получены изображения 
микроструктуры и поверхностного рельефа образцов после охлаждения в магнитном поле. Проанализированы 
кинетические изменения и зависимости объемных скоростей превращения от времени изотермической 
выдержки. Установлено, что действие постоянного магнитного поля напряженностью 1,6 МА/м увеличивает 
объемную скорость превращения в 1,808 раза (для стали 65Г) ив 1,687 раза (для стали 45Х). 

Обсуждение и заключение. Зафиксированы результаты наблюдений за изменением поверхностного рельефа 
при охлаждении без магнитного поля и в магнитных полях различной напряженности. Это позволило сделать 
вывод о стимуляции внешним магнитным полем бейнитного превращения вместо исходного перлитного. 
Микроструктурные изменения объясняются влиянием поля на магнитное состояние исходной фазы. 


Ключевые слова: магнитное стимулирование бейнитного превращения, перлитное превращение, улучшение 
свойств сталей, микроструктурные изменения сталей, вакуумное травление, магнитострикционные деформации 
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Введение. Как известно, магнитное поле влияет на термодинамику [1], в частности на равновесную 
температуру фазовых переходов [2]. Значительного смещения критических температур можно достичь лишь 
при использовании достаточно сильных (импульсных) магнитных полей [3], что сопряжено с определенными 
техническими трудностями при реализации в производственных условиях. Опыты свидетельствуют о 
появлении а-фазы в конструкционных и инструментальных сталях под действием полей напряженностью 
1,6-2,4 МА/м [1]. Эти процессы проходят при температурах, которые намного превышают равновесную и не 
укладываются в предполагаемый по термодинамическим оценкам сдвиг [4]. Влияние внешнего магнитного 
поля при термической обработке стали выражается в увеличении темпа [5] и дисперсности продуктов 
превращения [6]. Из [7] известно, что при промежуточном (бейнитном) превращении под влиянием магнитного 
поля напряженностью до 2 МА/м структурные изменения позволяют получать больший запас пластичности при 
сохранении или некотором повышении прочностных показателей. 

Имеющиеся в у-фазе области с ближним магнитным порядком при действии внешнего магнитного поля 
испытывают магнитострикцию. Это происходит из-за магнитной неоднородности аустенита [8] и вызывает 
изменение в поле упругих сил решетки, причем энергия образования критического зародыша @а-фазы 
уменьшается. Внешнее магнитное поле увеличивает количество и размеры ферромагнитных кластеров [1]. В 
результате при охлаждении в магнитном поле мультипликативно растет число центров зарождения. 

В [9] показано, как ближний магнитный порядок в аустените влияет на смену сценария фазового 
превращения. Из [10] известно, что магнитное состояние у-фазы определяет превращение в тот или иной 
продукт ферромагнитной а-фазы. Это может быть феррит, перлит, бейнит или мартенсит. 
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Цель представленной работы — экспериментальная проверка возможности смены перлитного сценария 
превращения на бейнитный при воздействии внешним магнитным полем до 1 МА/м. Кроме того, следует 
оценить изменение кинетики бейнитного превращения при наложении поля. 

Материалы и методы. Использовались образцы одной плавки из стали 65Г. Их химический состав 
контролировался на оптико-эмиссионном спектрометре МазеПап О8 (таблица 1). 

Таблица 1 


Среднее содержание элементов в образцах 
Массовая доля, % 
еря м | @ | 
65Г 0,65 0,20 0,97 0,21 0,009 0,0012 0,08 0,13 


Марка стали 


Полированные образцы помещали в вакуумную камеру установки для высокотемпературных исследований 
«ИМАШ 20-75». Данное оборудование обеспечивает разрежение в рабочей камере 1,3-10 - 6,6-10-° Па. Это 
позволяет реализовать метод вакуумного травления. Структура выявляется в результате испарения в вакууме 
при высокой температуре под влиянием поверхностного натяжения. При охлаждении в процессе фазовых 
превращений, сопровождающихся сдвиговыми процессами, на поверхности образца появляется 
соответствующий рельеф. 

Образцы подвергали резистивному нагреву до - 1000 °С в течение 10 минут и охлаждали с отводом тепла в 
водоохлаждаемые медные электроконтакты со скоростью -28-32 °С/с. В соответствии с диаграммой 
изотермического распада аустенита стали 65Г получаемая скорость охлаждения соответствовала пересечению 
«носа» области начала перлитного превращения. Высокая температура нагрева способствовала росту 
аустенитного зерна и вакуумному травлению его границ, что требовалось для видеозаписи поверхностного 
изменения рельефа в процессе фазового превращения. 

Для измерения температуры спай термопары приваривали в середине образца. Во время охлаждения 
протекающие на поверхности процессы регистрировали с помощью цифровой окулярной насадки ЕаК!пв. 
Полученные данные о микроструктуре поверхностного рельефа и внутри образца обрабатывали для 
определения объемной доли структурных составляющих. С этой целью задействовали аналитическое 
программное обеспечение ЗАМ 800 на 16-25 полях зрения. 

Внешнее магнитное поле создавали электромагнитом, интегрированным в вакуумную камеру 
установки [11]. В экспериментах задействовали поля напряженностью 400 кА/м и 800 кА/м. Образцы без поля 
обрабатывали тем же способом, но при извлеченном из вакуумной камеры электромагните. 

Результаты исследования. При бейнитном или мартенситном сдвиговом характере превращения на 
поверхности отполированного образца должен появляться рельеф. Если этого не происходит, речь идет о 
перлитном превращении. В условиях эксперимента скорость охлаждения была недостаточна для закалки на 
мартенсит стали 65Г. Можно было ожидать конкуренции между перлитными и бейнитным превращением в 
зависимости от наличия или отсутствия внешнего магнитного поля в процессе охлаждения. 

На рис.1 представлены скрины видеокадров изменения поверхностного рельефа, наблюдаемого при 
превращении переохлажденного аустенита в различные моменты времени (т). Каждой строке на рис. 1 
соответствует указанная напряженность внешнего магнитного поля (Н). 


Н, кА/м 


400 


800 


1875 2500 ме 


Рис. 1. Изменение поверхностного рельефа во время охлаждения в зависимости от напряженности внешнего магнитного 
поля (Н, кА/м) и времени превращения (т, мс) 
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На рис. 2 показаны микроструктуры поверхностного рельефа и внутри образцов после различных режимов 
обработки. 


Рис. 2. Микроструктура поверхностного рельефа (а, б, в) и внутри образца (г, д, е) после обработки: а, г — без поля; 0, д — в 
магнитном поле напряженностью 400 кА/м; д, е — в магнитном поле напряженностью 800 кА/м 


При обработке без поля изменения поверхностного рельефа наблюдались лишь в отдельных мелких 
образованиях (1-я строка, рис. 1 и рис. 2 а). На видеозаписи охлаждения зерен аустенита было чуть заметно 
распространение волнового процесса по поверхности. По-видимому, это отражение фазового превращения. 
Микроструктурный анализ внутренних слоев образца показал наличие подавляющего количества (98%) 
перлитной структуры (рис. 2 г). 

С наложением внешнего магнитного поля поверхностный рельеф формировали зарождающиеся и растущие 
бейнитные иглы (2-я строка, рис. 1 и рис. 2 6). Эти кристаллы нельзя отнести к мартенситным из-за их 
медленного роста (рис. 1). С увеличением напряженности поля повышалась интенсивность сдвигового 
превращения на поверхности (3-я строка, рис. 1 и рис. 2 в), образовывался пакетный бейнит. 

При действии магнитного поля напряженностью 400 кА/м объемная доля бейнита увеличилась 
втрое — до 6 % (рис. 2 д) по сравнению с обработкой без поля. В магнитном поле напряженностью 800 кАм 
объемная доля бейнита составляла уже 8-10 %. Это позволяет сделать вывод о стимулировании смены 
сценария превращения из перлитного в бейнитное при обработке во внешнем магнитном поле. 

Если при охлаждении исходную фазу (аустенит) сохранять как можно дольше, то усиливающийся в ней с 
понижением температуры ближний магнитный порядок (магнитная неоднородность) в своей критической 
степени приведет к атермическому мартенситному сценарию превращения. Это известно из [9]. В [1] показано, 
что наложение внешнего магнитного поля в процессе охлаждения аустенита приводит к дополнительному, 
вынужденному магнитному расслоению аустенита из-за увеличения количества, размеров и времени 
устойчивого существования ферромагнитных кластеров. В условиях текущего эксперимента искусственно 
усиленный внешним полем ближний магнитный порядок в у-фазе способствует протеканию бейнитного 
превращения вместо перлитного, естественного для данных условий охлаждения без наложения поля. 

Под действием магнитного поля значительно меняется кинетика бейнитного превращения. Об этом 
свидетельствует отмеченное в ходе эксперимента усиление бейнитной реакции при увеличении напряженности 
магнитного поля. Кроме того, в [7] рассматривалось наложение внешнего магнитного поля напряженностью 
1,6 МА/м при изотермической выдержке различных сталей. Итоги указанной работы также подтверждают 
сделанное выше заявление о кинетике превращения. 

Для выяснения механизма влияния внешнего магнитного поля на кинетику промежуточного превращения 
учитываются: 

— специфика роста бейнитных кристаллов, который зависит от скорости отвода атомов углерода от у/о границы; 

— возникающие при превращении структурные напряжения, обусловленные изменением удельного объема 
превращающихся фаз. 
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Данный подход обусловлен тем, что процессы релаксации у межфазной границы сильно заторможены при 
пониженных температурах [12]. В такой ситуации градиент напряжений вызывает дрейф атомов углерода 
вопреки концентрационной неоднородности. Это играет определяющую роль в скорости роста нижнего 
бейнита. Перемещение атомов С ведет к снижению его концентрации в объемах у фазы по фронту роста 
пластин бейнита. В результате возникает контролируемый концентрационным градиентом диффузионный 
поток. Он направлен навстречу растущему кристаллу и уменьшает скорость его роста. Зависимости, 
описывающие эти процессы, есть в работах Л.Н. Александрова и Б.Я. Любова [13]. 

Под действием внешнего магнитного поля возникают магнитострикционные деформации. Их упругая 
энергия может вносить определенный вклад в энергию взаимодействия диффундирующих атомов с полем 
структурных напряжений. Следует ожидать усиления дрейфа углерода и увеличения скорости роста кристаллов 
а фазы. Скорость роста по механизму дрейфа в соответствии с расчетами Л.Н. Александрова и Б.Я. Любова 
описывается зависимостью: 


и 2РБ , (1) 


вт 
ВАТ 1-| —®. 
Со 


Здесь р — коэффициент диффузии углерода в аустените; Ак, — радиус критического зародыша при данном 
переохлаждении; Снъ и Со — соответственно концентрация углерода в а фазе и исходном аустените. Р 
характеризует энергию взаимодействия диффундирующих атомов с полем структурных напряжений, 
вызванных дилатацией = при изменении объема в процессе превращения, и находится из соотношения: 


8 
Р=—л’® . 2 
9 (2) 


где г. — радиус атома углерода; © — параметр, характеризующий в твердых растворах внедрения зависимость 
постоянной решетки от концентрации атомов растворенного элемента; Ё — модуль Юнга; у — коэффициент 
Пуассона; и = (Ё/2) : (1-у) — коэффициент Ламэ (модуль сдвига по направлению [100]); т = (а-Б) - (а+), 
гдеаи Ь — размеры полуосей кристалла в форме эллипсоида вращения. 

Расчет по формуле (1) проводили для углеродистой стали при: Си= 0,7; т = 0,9; у=0,3; и = 73,5 Гпа; 


—131 2 
р о моль) г ; = 0,01; © = 0,02. Величину Акр определяли из соотношения 


Кр =26 / АРо. Поверхностное натяжение на когерентной границе с = 0,2 Дж/м?; АЕо при Т= 600 К составляет 
315 МДж/мз. Значение С.„ъ=0,4 определяли как концентрацию, соответствующую температуре начала 


Р=0,0999 с 
с 


мартенситного превращения (600 К). Найденная таким образом скорость роста кристаллов а фазы составляет 
— 8-10-6 см/, т. е. является величиной того же порядка, что и Их» -- 10-6 см/с, полученная экспериментально для 
стали 65Г при видеосъемке изотермического образования рельефа. 

При включении магнитного поля измерение скорости роста дало те же результаты, т. е. экспериментально 
не обнаруживается влияние магнитострикционных напряжений на интенсификацию дрейфа. 

Допустим, упругое смещение равно истинной магнитострикции парапроцесса /, = 0,5:10-4 [14]. В этом случае 
при Н = 1,6 МА/м и [= 400 °С напряжения от магнитострикции составляют о), = Е› = 10 Мпа. Они вносят весьма 
малый вклад в энергию взаимодействия диффундирующих атомов с полем структурных напряжений. Это 
объясняется тем, что величина магнитострикции на два три порядка меньше величины дилатации 
(0,01 = 0,07 [15]) при сдвиговом у-—а переходе. По этой причине оценка вклада магнитострикции по 
формуле (1) дает исчезающе малую разницу в значениях скорости роста при обработке без поля и в магнитном поле. 

Ускорение процесса образования бейнита можно фактически оценить при помощи уравнения 
А.Н. Колмогорова [16], устанавливающего связь между степенью превращения, скоростью зарождения и 
ростом центров новой фазы. Изменение объемной скорости превращения в магнитном поле из-за уменьшения 
работы образования ферромагнитного зародыша критического размера [1] можно найти по формуле: 

он "| 


= ехр 1+ - г 
уу) | уй 


(3) 


где о — объемная скорость превращения (индексы 0 и Н соответственно обозначают условия обработки без поля 
и во внешнем магнитном поле); И’ — энергия образования равновесного зародыша; / — намагниченность; 
Н — напряженность магнитного поля; АХ — удельная «химическая» движущая сила; Д/” — свободная энергия 
образования одного ферромагнитного кластера [1]. 
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Согласно классической теории Л. Кауфмана и М. Коэна [17], энергия зарождения центров новых фаз равна: 


27.0202 
и _ 95 ра (4) 
(м+м) 
Здесь ДГ — работа образования ферромагнитноупорядоченных кластеров в аустените; 9 — параметр, 


учитывающий влияние энергии упругих деформаций; о — поверхностное натяжение. Из [17] известно, что 
0203 = 9,92 . 10 Джз/м®. 

Вычисления проводились для температур 543 К (сталь 65Г) и 628К (сталь 45Х) при напряженности 
магнитного поля Н = 1,6 МА/м. Энергия образования равновесного зародыша И’ определялась при удельной 
химической движущей силе Д/= 150 МДж/м? с учетом напряженности поля и среднего размера ферромагнитного 
кластера - 1,8 нм [1], для которого получена величина ДХ = 0,63 МДж/м3. Расчет по формуле (3) дает для 
бейнитного превращения в стали 65Г он / го = 1,804, а в стали 45Х он / и = 1,665. 

Зависимости объемной скорости превращения (рис.3) получены с использованием экспериментальных 
данных [7] о степени бейнитного превращения в сталях 65Г и 45Х при различных временах изотермической 
выдержки без поля и в магнитном поле напряженностью 1,6 МА/м. 

1,5 3,00 


ю 
[1] 
‘л 


1,0 


1,5 


Объёмная скорость 
превращения, %/с 


Объёмная скорость 
превращения, %/с 


0 100 200 300 0 50 100 150 200 
Время, с Время, с 


а) 6) 


Рис. 3. Зависимости объемной скорости превращения (вертикальная ось, %/с) от времени изотермической выдержки: 
а — для стали 65Г при температуре 375 °С; 6 — для стали 45Х при температуре 320 °С. Синяя кривая отражает данные по 
обработке без поля, красная — в магнитном поле. На каждой линии отмечены экспериментальные точки 


Рис. 3 позволяет рассмотреть отношение объемной скорости превращения при обработке в поле к обработке 
без поля он/о в самом начале бейнитной реакции (экстраполированное значение в т=1 с). Для стали 65Г 
получится 1,808, для стали 45Х — 1,687. Это близко к приведенным выше теоретическим оценкам: 1,804 для 
стали 65Г и 1,665 для стали 45Х. В процессе превращения скорости меняются. Если рассмотреть самые первые 
экспериментальные точки на графиках, то: 

— для стали 65Г при т= 40 сон/е = 1,75; 

— для стали 45Х при т= 28 сон/ш = 1,629. 

При более длительной изотермической выдержке отношения объемных скоростей уже мало коррелируют с 
расчетными. 

Заключение. Экспериментально установлено, что, задействуя внешнее магнитное поле, можно сменить 
сценарий превращения с перлитного на бейнитное. Это объясняется усилением степени ближнего магнитного 
порядка в аустените в результате увеличения количества и размеров ферромагнитных кластеров в у фазе под 
влиянием магнитного поля. С ростом его напряжения увеличивается магнитная гетерогенность аустенита. Как 
следствие, повышается скорость бейнитного превращения. Расчетные значения увеличения объемной скорости 
превращения при действии внешнего магнитного поля хорошо согласуются с экспериментальными данными 
для стартовых скоростей бейнитного превращения. 
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